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Stresses in Laminated Composites 

ABSTRACT 

An approximate solution for interlaminar stresses in finite width, laminated composites sub- 
jected to uniform extcnsional and bending loads is presented. The solution is based upon the 
principle of minimum complementary energy and an assumed, statically admissible stress state, 
derived by considering local material mismatch effects and global equilibrium requirements. The 
stresses in each layer are approximated by polynomial functions of the thickness coordinate, mul- 
tiplied by combinations of exponential functions of the in-plane coordinate, expressed in terms of 
fourteen unknown decay parameters. Imposing the stationary condition of the laminate comple- 
mentary energy with respect to the unknown variables yields a system of fourteen non-linear 
algebraic equations for the parameters. Newton’s method is implemented to solve this system. 
Once the parameters are known, the stresses can be easily determined at any point in the lam- 
inate. 

Results are presented for through-thickness and interlaminar stress distributions for angle-ply, 
cross-ply (symmetric and unsymmetric laminates), and quasi-isotropic laminates subjected to uni- 
form extension and bending. It is shown that the solution compares well with existing finite ele- 
ment solutions and represents an improved approximate solution for interlaminar stresses, pri- 
marily at interfaces where global equilibrium is satisfied by the in-plane stresses, but large local 
mismatch in properties requires the presence of interlaminar stresses. Further, the contributions 
of both global equilibrium and local material mismatch effects to the stress field are clearly del- 
ineated. The results indicate that the significance of local mismatch effects is dependent on lam- 




inate stacking sequence. The demonstrated accuracy and efficiency of the solution make it 
ideally suited for parametric studies. 
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CHAPTER 1 

INTRODUCTION 

1.1 Introductory Remarks 

The high strength-to-weight, and stiffness-to-weight ratios of composite materials and their 
tailorability to meet strength and stiffness requirements has led to the increased use of 
composites for structural designs, particularly in aerospace applications. With this increased use 
has come significant interest in the failure mechanisms of composite materials. Because 
composite materials are heterogeneous and anisotropic, failure modes occur that are quite 
different from those seen in more conventional isotropic materials. Laminated composites 
exhibit two basic failure modes: 1) in-plane fracture, and 2) out of plane delamination * ! * failure, 
matrix failure, fiber/matrix debonding or fiber splitting 121 

Experimental studies have shown that the mode of failure and ultimate failure load of 
laminates are dependent upon the laminate stacking sequence and the layer thicknesses.* 2 ' 8 * 
This phenomenon cannot be explained by classical laminated plate theory (which predicts a 
planar stress state) combined with in-plane fracture theories, but is attributed to the presence of 
interlaminar stresses near the free edges of composite laminates. 

Interlaminar stresses arc caused by the mismatch, or difference, in the material properties 
between the individual laminae and the laminate and the mismatch in properties between 
adjacent laminae in the presence of a free edge. Individual layers of a laminate will deform 
differently, when subjected to the same axial strain, because of differences in their material 
properties. In a laminate, however, the layers are bonded together, and displacement continuity 
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at the layer interfaces requires development of interlaminar stresses to equalize the differential 
deformations and to maintain equilibrium. A detailed discussion of the mechanics of free edge 
stresses is provided in reference [9]. 

The interlaminar stresses a 33 , a 23 , and (j) 3 . shown in Figure 1.1, act upon planes parallel to 
the interfacial planes between laminae. They exist only within a very local region near the free 
edges of a laminate and are therefore known as a boundary layer effect or free edge effect. For 
fiber reinforced composites the interlaminar stresses are transferred between plies through the 
matrix material that bonds them together. This interfacial region is relatively weak and if the 
interlaminar stresses are high enough the laminated structure will fail, due to delamination, at 
loads much lower than those predicted by in-plane failure theories. 

A necessary tool to aid in understanding and ultimately preventing delamination type failures 
is an efficient analytical method which provides reasonably accurate stress predictions in the 
boundary layer region. The need for such a method is particularly acute in design stages to 
avoid delamination prone laminates when a large number of possible structural configurations 
have to be evaluated quickly and economically. Numerous investigators have proposed a variety 
of methods for calculating interlaminar stresses. The majority of these solutions are numerical in 
nature and are plagued by computational limitations, particularly with regard to memory 
requirements. Consequently they become intractable when the number of layers in a laminate 
becomes even moderately large or when calculations have to be performed repetitively as in an 
optimization process. 

In practical applications composite panels may consist of many layers (100 layers in aircraft 
structures is not uncommon) of different orientations, thicknesses and material properties. Thus, 
design of even the simplest composite structural component may involve a large number of 
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design variables in addition to a large number of design constraints. These factors make 
composite panel design an ideal candidate for numerical optimization. 

Several papers, a few of which are listed in the references.^ 10 ' 18 ' have been published on the 
use of numerical optimization for designing composite structures. To the author’s knowledge, 
however, none of the published work on composite design using mathematical optimization 
techniques considers interlaminar stresses in the problem formulation. This gap in the literature 
apparently exists because of the complexity and computational inefficiency of the majority of the 
methods currently available for predicting interlaminar stresses. 

1.2 Objective and Scope 

The above discussion suggests the need for more efficient approaches for calculating the 
three dimensional stress field near free edges in laminated composites. An analytical approach is 
preferred so that the method can be incorporated into a design process which uses numerical 
optimization techniques. The objective of this research is to provide an approximate analytical 
model for laminate stress analysis and demonstrate its usefulness. The approximate model is 
developed for a finite width symmetrically or unsymmetrically laminated coupon with straight 
free edges subject to uniform extensional or bending loads as shown in Figure 1.1. These loads 
are considered because they are common in practice. Combined loads can then be analyzed by 
superposition. The straight, free edge coupon was chosen because it is the simplest 
configuration to analyze and there are numerous results available in the literature that can be 
used to verify the model developed. Laminates with many plies as well as hybrid laminates can 
be analyzed. Also, although the formulation that follows is presented with reference to the 
simple plate shown in Figure 1.1, more complicated structural configurations can be analyzed 
using the methodology developed provided the in-plane stress field in the interior can be 
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obtained from an analytical solution or from some general analysis technique. The methodology 
can then be employed in a global-local analysis to obtain refined stress solutions in regions of 
high stress gradients with a coarser global solution used to define the response outside of these 
regions. 

The remainder of this thesis is divided into several sections. Chapter 2 includes a literature 
review of various methods for predicting interlaminar stresses and a brief discussion of the 
present method. In Chapter 3 the analytical model is developed. The approach is an extension 
of previous work by Kassapoglou and Lagace I 1 - 19 " 21 ! and is based upon an assumed stress state 
and the principle of minimum complementary energy. The methodology reduces the stress field 
determination for a general laminate to the simultaneous solution of 14 non-linear equations. 
Newton’s method is implemented to solve this system. In Chapter 4 the solution is verified and 
its advantages and limitations are identified by comparison with existing finite element and 
analytical model results. Comparisons are made for angle-ply, cross-ply and more general 
laminate configurations. In Chapter 5 additional results are presented demonstrating the utility 
of the technique and the effect of load conditions and stacking sequence on interlaminar stresses. 
Finally, Chapter 6 closes with a summary and conclusions of the study. Recommendations for 


future work are also provided. 
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CHAPTER 2 


LITERATURE REVIEW 


Interlaminar stresses have been studied for over twenty years, in hundreds of articles. Only a 
few are discussed here. For a more complete coverage see the review article by Salamon . 1221 The 
majority of the work has concentrated on the analysis of symmetric laminates with straight free 
edges subjected to uniform axial extension. A brief review of the more significant contributions 
to the understanding of interlaminar stresses in these laminates is presented first. This is 
followed by a review of the literature pertaining to interlaminar stress calculations for laminates 
in uniform bending. The chapter closes with a summary and discussion of the present solution 
methodology. 

2.1 Laminates in Uniform Axial Extension 

The earliest investigations of interlaminar stresses were performed by Hay&shi , 1231 Hayashi 
and Sando 1241 and Puppo and EvcnsenJ 251 All of these researchers modeled the laminate as a set 
of anisotropic layers separated by isotropic shear layers. Their analyses neglected the 
interlaminar normal stress component O 33 and predicted a sharp rise in the interlaminar shear 
stress a 2 3 at the intersection of an interface and the free edge. 

In the same year, Pipes and Pagano performed the first numerical study of edge stresses in 
composite laminates . 1261 They studied the elastic response of a [+45/— 45], laminate subjected to 
a uniform axial extension (Figure 1.1). Noting St. Vcnant’s principle, they assumed the stresses 
to be independent of the axial coordinate Xj in regions away from the areas of load 
introduction. Under this assumption the general form of the displacement field is 
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u (Xj X 2> X 3 ) = X,e„ + U(X 2 X 3 ) 

(2.1) 

v(X,X 2> X 3 ) = V(X 2i X 3 ) 

(2.2) 

w(X, X 2 X 3 ) = W(X 2 X 3 ) 

(2.3) 


where u, v, and w are displacements in the X x , X 2 , and X 3 coordinate directions respectively, 
and en is the applied axial strain. The reduced elasticity equations governing the laminate 
behavior were then formulated and solved using the finite difference method. 

Their results showed a planar stress field over most of the laminate in agreement with 
classical lamination theory (CLT). In regions near the laminate free-edge, the lamination theory 
results were perturbed by the presence of the interlaminar stress components, a 13 , a 23 , and a 33 . 
The interlaminar stresses were shown to decay rapidly with distance from the free edge and were 
zero outside of a region of width approximately equal to the laminate thickness. Therefore they 
concluded that interlaminar stresses are a boundary layer or an edge effect. 

Three of the predicted stress components, O23, cr 22 and a 33 were very small while the 
interlaminar shear stress Oi 3 was quite large. In addition they noted that the magnitude of Oj 3 at 
the intersection of the interface and the free edge increased with increasing grid refinement. 
Based upon these results and those of Bogy^ 271 and Hess^ 28,29 1 for bonded quarter planes of 
dissimilar materials, they concluded that Oi 3 is singular at this point. The results of this model 
along with some simplified models for predicting selected interlaminar stress components were 
used to explain the relationship between interlaminar stresses and the differences in 
experimentally observed strengths of similar laminates.^ 4 ’ 6,30 ^ 

Finite element solutions soon followed. The first finite element solution was provided by 
Isakson and Levy. 131 ' They used a displacement based formulation and like Puppo and Evcnsen 
modeled the laminate as a combination of anisotropic layers separated by isotropic shear layers, 
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thus neglecting the interlaminar normal stress. They analyzed a [445/— 45], laminate and their 
predictions for the interlaminar shear stress O 13 agreed well with those of Pipes and Pagano. 
Rybicki 1321 used a three-dimensional finite element analysis to obtain approximate solutions for a 
symmetric laminate subjected to in-plane loading. He used a complementary energy formulation 
with assumed stress states derived from the three-dimensional Maxwell stress functions. His 
formulations provided predictions for all three interlaminar stresses and the results showed good 
agreement with those of Pipes and Pagano. 

Later, to improve solution efficiency, several investigators adopted the Pipes-Pagano 
approach and solved the tensile coupon problem using quasi-three dimensional formulations. 
The first two-dimensional finite element analysis for the quasi-three dimensional problem was 
conducted by Herakovich et. al to study mechanical and thermal edge effects in cross-ply and 
angle-ply laminates .^ 33,341 Application of this type of formulation to additional laminate 
configurations soon followed .^ 35 ' 391 Wang and Crossman 1351 analyzed 5 laminate configurations; 
two cross-ply laminates, an angle-ply laminate, and two quasi-isotropic laminates. By invoking 
a skyline storage scheme, they were able to use a much finer mesh than had been used in 
previous analyses and obtained a more accurate description of the stress field in the vicinity of 
ply interfaces and the free edge. They noted that the interlaminar normal stress 033 is also 
singular at these points. 

The finite element method provided a means for obtaining solutions for a variety of laminate 
configurations and geometries. Numerous solutions were obtained. These solutions greatly 
increased the understanding of the free edge problem and the mechanisms contributing to 
interlaminar stress development. The limitations of numerical procedures for laminate stress 
analysis, however, also became evident. First, it was quickly realized that the numerical 
solutions are not economical. Because of the singular nature of the problem, extremely fine 
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meshes or fine finite difference grids are required in edge regions in order to obtain reasonably 
accurate predictions of the field variables. Wang and Crossman, 1351 for instance, used 192 
elements (16 through the thickness) to model each ply in a four ply laminate. Meshing 
requirements like these and the resulting computer memory and time requirements made 
analysis of practical laminates prohibitive. Second, it was found that results very near the free 
edge obtained using different formulations were not consistent. Different researchers not only 
predicted different magnitudes for interlaminar stress components in this region but also, in 
some cases, predicted different signs. Some attributed this anomaly to improper satisfaction of 
the free edge boundary conditions. 1371 Others suggested that lack of symmetry in the stress tensor 
at the singularity 1401 may be the cause of the inconsistencies. 1381 

In an effort to resolve these inconsistencies and develop more efficient reliable methods for 
laminate stress analysis, several analytical solution methods were proposed. One of the first 
fairly sophisticated approaches which was capable of predicting both interlaminar shear stress 
components as well as the interlaminar normal stress component was provided by Tang 1411 and 
Tang and Levy. 1421 They extended the boundary layer theory of plane stress of isotropic elasticity 
developed by Reiss and Locke 1431 to the analysis of laminated composites. Using a zeroth order 
approximation to the boundary layer problem, the solution in the boundary layer region was 
separated into a torsion problem and a plane strain problem. The solution exactly satisfied the 
equilibrium equations and compatibility equations but some of the boundary conditions were not 
satisfied or were satisfied in an average sense. A similar approach was used by Hsu and 
Herakovich 1441 in the study of angle-ply laminates. Using a perturbation method, they matched 
an interior solution, where classical lamination plate theory is assumed to be valid, to a boundary 
layer solution. Their results suggested that both interlaminar shear stress components Oi 3 and 
a 2 3 arc singular in the boundary layer region. Another approach was presented by Wang and 
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Dickson.* 451 They expressed displacements and the interlaminar stresses in each layer in a scries 
of Legendre polynomials, and used Galcrkin’s method to obtain a system of equations for the 
unknown constants. Out of plane warping was neglected in the displacement assumptions so 
their model was only applicable to cross-ply laminates. They stated that the method is capable 
of handling laminates with a large number of plies but didn’t provide any results to support this 
claim. Also because of convergence difficulties (for larger b/t ratios) in stresses at an interface 
and the free edge, the solution is limited to very thin laminates. 

A similar approach was proposed by Bar-Yoseph and Pian. 1461 In their solution, the edge 
layer stress field is constructed using Legendre Polynomials that exactly satisfy the equilibrium 
conditions, traction continuity conditions, and the stress free edge conditions. Later they 
incorporated their assumed stress states into a mixed hybrid finite element formulation. 1471 This 
method was extended by Bar-Yoseph and Siton to include nonlinear material behavior. 1481 

Another variational approach was provided by Pag ano. 1491 He proposed an approximate 
solution based upon the extension of Reissner’s theorem 1501 to a laminated body. Requirements 
for an acceptable laminate field theory were established; all stress components are non-zero, 
displacement and traction continuity are satisfied at all interfaces, and each layer or sublayer 
(more than one sublayer per ply is permitted) is in equilibrium. In establishing layer 
equilibrium, free edge conditions are imposed on force and moment resultants rather than on 
point-wise tractions. A model which is based upon assumed stress fields in each layer is 
developed which satisfies these criteria. Explicit functions are assumed for the through thickness 
variations of stresses. Minimization of Reissner’s functional over the entire laminate results in a 
system of 13N differential equations, where N is the number of sublayers in the laminate, for the 
in-plane variations of the field variables. 
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Pagano 1511 also delineated the theory for the analysis of an axial coupon. The field equations 
in this case arc constant coefficient, linear differential equations, in the width coordinate y. The 
homogeneous solution for each independent variable is then a sum of exponential terms of the 
form 

f = Fe Xy (2.4) 

The method provides accurate stress distributions in regions near free edges but very large 

magnitudes of X, obtained for large N, limited solutions to laminates with N < 6 because of 

computer overflow/underflow violations. Also, although the solution does not include a 

singularity, the o 13 and o 33 stress components increased as the number of sublayers used per 

layer increased. This behavior is similar to that observed with increased mesh refinement when 

using the finite element method. 

Pagano and Soni [52 ’ later took this model and, using a global-local variational formulation, 
developed a ply/sub-laminate analysis. In regions where a detailed response is required (local 
region) each ply is represented by the model described above. The remaining areas, e.g. sub- 
laminates or global regions, are represented by effective elastic properties. 1531 The method shows 
promise but appears to be somewhat sensitive in its predictions to the choice of the global and 
local domains. Also, if stresses are desired at each interface of an N layered laminate, the 
global/local analysis must be exercised several times with the local domain containing the 
interface of interest. Rehfield et. al 154,551 employed a similar approach using their refined 
theories for the behavior of anisotropic plates for the ply/sublaminate models. 156,571 Their 
method results in a set of 8N-3 equations, where again N is the number of sublayers. 

Most of the solutions described above suggest the presence of a singularity at the intersection 
of an interface and the free edge. To incorporate the singularity in the formulation, the nature of 
the singularity must be known before hand. Up to this point Bogy’s 1271 work on the singularities 
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on isotropic quarter planes had not been extended to anisotropic materials. Consequently, to 
better understand the boundary-layer effect in composite laminates, some investigators set out to 
determine the exactly the singularity at the intersection of the interface of two laminae and the 
free edge. One such study was presented by Wang and ChoiJ 58,59 ^ Their formulation is based 
on the theory of elasticity and Lekhnitskii’s^ 60 ' complex stress potentials and leads to a pair of 
coupled governing partial differential equations. The homogeneous solution to the equations is 
obtained using an eigenfunction expansion. The homogeneous solution showed the existence of 
a singularity of the form y -5 at the intersection of an interface and the free edge. They found 
that the order of the singularity 5 is in general very weak and is dependent only on the material 
constants and fiber orientations of plies adjacent to the interface of interest. Similar studies by 
Zwiers, Ting and Spilkcr, 16 ' 1 and Dempsey and Sinclair^ 62,63 ! showed that singularities of the 
form ln(y), On(y)) 2 , (ln(y)) 3 , etc. are also present for some combinations of adjacent layers, in 
addition to the y singularity. These results, along with numerical studies on the 
singularities, [391 arc significant because they showed that although mathematical singularities 
exist, they are generally very weak and act over such small distances that approximate solutions 
that do not incorporate the singularity are accurate except in regions very near the free edge. 

2.2 Laminates in Uniform Bending 

Few studies have been conducted on laminates in bending. Salamon [641 presented a solution 
for finite width laminates uniformly bent by end moments applied about the X 2 -axis (see Figure 
1.1). Using an approach similar to that of Pipes and Pagano^ the elasticity equations are 
formulated and solved using the finite difference method. He finds that the interlaminar shear 
and normal stress distributions are similar for a laminate in uniform bending to those of a 


laminate in uniform extension. 



13 


Finite element studies were conducted by Murthy and Chamis^ 65 ' for a variety of load 
conditions including in-plane and out-of-plane bending, and by Chan and Ochoa^ 66,67 ^ for 
laminates under torsion and bending loads. Kassapoglou* 68 * extended his analysis for extension 
loading* 20 * to combined loading cases and bending. 

2.3 Summary and Discussion 

The above discussion gives an indication of the variety of solutions proposed for obtaining 
free edge stresses in laminated composites. These solutions have increased the understanding of 
the mechanics of interlaminar stress development and their effects on the performance of 
laminated composite structures. Most of these solutions, however, are constrained 
computationally by the size of laminate system they can handle and therefore have limited 
practical application. Apparently some tradeoff needs to be made between solution accuracy of 
the complicated solutions described above and solution efficiency offered by simplified 
approximate models. 

Recently Kassapoglou and Lagace* 21 * proposed a simplified, approximate technique for 
determining the stress field in the vicinity of straight free edges of a laminated coupon. A very 
similar approach was presented by Engrand . 1691 Kassapoglou and Lagaces’ analysis is based on 
the principle of minimum complementary energy and an assumed stress state obtained by 
considering global equilibrium requirements. Generic forms of stress distributions that exactly 
satisfy the equations of equilibrium, the traction continuity conditions and the free edge stress 
boundary conditions are assumed within each layer. Explicit polynomial expressions arc used 
for the through-thickness variations of the stresses, while the in-plane variations are taken to be 
combinations of two decaying exponential functions expressed in terms of two unknown decay 
parameters. The unknown decay parameters are functions of ply material properties, orientation 
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and thickness, and laminate stacking sequence, and arc determined by minimizing the laminate 
complementary energy. Once these parameters are determined the stresses can be evaluated at 
any point in the laminate. 

A distinguishing feature of Kassapoglou and Lagaces’ model is that run times and computer 
memory requirements are a linear function of the number of layers in the laminate. This makes 
for an extremely efficient design tool that can be used to analyze laminates with many layers. 
The model, however, has trouble predicting interlaminar normal stresses and in some cases 
interlaminar shear stresses. 

Kassapoglou and Lagaces’ (KL) method serves as a basis for the approximate solution 
developed in this investigation. In their stress assumptions only the mismatch between laminae 
and laminate material properties is considered. The improved solution includes additional terms 
in the stress assumptions which account for the effect of mismatch in engineering properties 
between adjacent layers of a laminate. Specifically, the mismatches in coefficient of mutual 


influence and Poisson’s ratio are considered. 
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CHAPTER 3 


SOLUTION FORMULATION 


The problem considered is the uniform axial extension or bending of a multi-layered 
laminated plate. Interest focuses on calculating the stress field in the vicinity of the free edge, 
i.e. in the boundary layer region. In this chapter, the mathematical boundary value problem for 
determining the free edge stresses is formulated from the linear theory of elasticity. Because of 
the analytical complexities of the three dimensional elasticity equations that must be solved 
within each layer, coupled with the requirement of continuous displacements and stresses at 
interfaces between layers, an exact elasticity solution for stress analysis in practical laminates is 
not feasible. 158,591 Hence, an approximate solution is proposed. The approximate solution is 
based upon the principle of minimum complementary potential energy and stress assumptions 
constructed in such a manner as to simplify the equations to be solved, while retaining the 
necessary three dimensional characteristics of the stress field. The stress assumptions exactly 
satisfy all of the equilibrium requirements. The compatibility equations and displacement 
continuity conditions are satisfied in an average sense through minimization of the laminate 
complementary energy. 

A singularity is not included in the stress assumptions. As previously mentioned, previous 
investigators 158,61,621 have shown that a very weak stress singularity is present near the 
intersection of interfaces and the free edge of composite laminates. However, as Pagano has 
noted, 1491 these singularities are artifacts of the effective modulus approach and do not exist in 
real materials. Also, the singularity is so weak and acts over such a small portion of the laminate 
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near the free edge that the assumption of material homogeneity on which the analysis is based 
breaks down. Further, Pagano suggested that when interpreting stress predictions using an 
effective modulus approach, average stresses, rather than point stresses, in regions of steep 
gradients may lead to more realistic conclusions regarding physical behavior. 170 ' These 
comments suggest that the stress singularities may be of only academic concern, and a solution 
that does not include a singularity is equally as valid as one that does, particularly in design 
applications where a qualitative comparison of the interlaminar stress severity in candidate 
laminates is the primary interest. 

The remainder of the chapter is devoted to development of the approximate solution. In 
Section 3.1, the elasticity problem is formulated. Section 3.2 summarizes classical lamination 
theory and the stresses of classical lamination theory are derived in terms of the mismatches in 
laminate and laminae material properties. These stresses and local mismatch considerations are 
then used to formulate the refined approximate model as discussed in Section 3.3. 

3.1 Problem Statement 

The geometry of a long, symmetrically or unsymmetrically laminated plate of finite width is 
shown in Figure 1.1. The laminate is built up of several layers reinforced by a system of parallel 
fibers oriented at an angle 0 with respect to the laminate longitudinal axis. Perfect bonding 
between adjacent layers is assumed. The laminate is assumed to be long enough so that away 
from the ends, where the loads are applied, the stresses and strains are independent of the axial 
coordinate. Another assumption made is that away from the edges the laminate is in a state of 
plane stress with the response defined by the classical lamination plate theory model. This 
assumption places a limitation on the geometry of laminates that can be accurately analyzed 
using the approximate solution and will be discussed further in Section 3.3.5. Also, 
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unsymmetric laminates will deflect out-of-plane when subjected to in-plane loads because of the 
membrane- flexural deformation coupling behavior they exhibit. These out of plane deflections 
are assumed small, so that geometric coupling effects can be ignored. 


An exploded view of the laminate, showing the in-plane and out-of-plane stress components 
is provided in Figure 3.1. The origin of the global coordinate system (X 1 ,X 2 ,X 3 ) is located at 
the center of the laminate, with the Xj, X 2 , and X 3 axes taken in the axial, transverse, and 
thickness directions respectively. Local coordinates x^\ and are established in each 
layer, where y (k) = b - x™ is measured from the free edge and z (k) is measured from the bottom 
of the kth ply. Beginning at the top surface of the laminate, the layers are numbered 
consecutively from 1 to N. The layers may have different thicknesses and may be different 
materials. Each layer is represented by a macroscopically homogeneous, linearly elastic, 
orthotropic material. Since the fiber axes of the individual layers are rotated through an angle 0 
with respect to the laminate axis, the material behavior of each laminae appears monoclinic in 
the global coordinate system. The constitutive equations for each layer then have the form 
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where [Sy ] is the transformed compliance matrix. 

3.2 Classical Lamination Theory 

3.2.1 Assumptions and Constitutive Relations 


According to classical lamination theory, the plate in Figure 1.1 acts as a single integral unit 
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with smeared elastic properties, and deforms under load in accordance with the Kirchoff 
deformation assumptions for thin plates. The elastic non-homogeneity of the laminate is taken 
into account by calculating the stresses in the individual layers using the laminate strains 
determined with these assumptions. A state of plane stress is assumed in each ply. The plane 
stress assumption implies 


= a g = a# = 0 


(3.2) 


and the stresses in the k* layer are given by 
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where Q,j are the reduced stiffnesses in the laminate coordinate system. These relations are 
used in conjunction with the Kirchoff deformation assumptions to define integrated laminate 
properties. From the Kirchoff assumptions for thin plates, the laminate strains are 


(e)x - {e°)x + X3{k)x (3.4) 

where {e°} x are the laminate middle surface strains and {tc} x are the laminate middle surface 
curvatures. Substituting the through-thickness strain variations (3.4) into the layer constitutive 
relations (3.3), yields expressions for the stresses in the k* 1 * layer in terms of the laminate middle 
surface strains and curvatures: 
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The laminate constitutive relations are then obtained by integrating equation (3.5) through the 
laminate thickness. This yields 
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where the laminate stiffnesses are 


f N l 
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+H 

[A,B,D]= J [Q] k [l,X 3 ,X|] dX 3 

-H 

and the force and moment resultants acting on the laminate are 
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(N) = J {d} k dX 3 

-H 
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(M) = J {dj k X 3 dX 3 
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(3.6) 


(3.7) 


(3.8a) 

(3.8b) 


3.2.2 CLT Stresses from Global Mismatch Considerations 


Stresses develop in the classical lamination theory because of the mismatch in material 
properties between the laminate and the individual layers comprising the laminate. This type of 
mismatch will be referred to as the global mismatch in material properties. To show the 
relationship between global mismatch and the classical lamination theory stress components 
consider a symmetric laminate subjected to in-plane uniaxial extension . For this loading and 
geometry all terms in the [B] matrix are zero, and the extensional response uncouples from the 
bending response. The laminate constitutive relations are then 

(N) = [AHe°) x (3.9) 


and the stresses in the kth ply are 
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where e°i, £22. and y?2 arc the laminate strains and arc constant throughout the laminate 
thickness. The laminate strains e-q and y?2 can be related to the applied strain e?i by the 



21 


laminate Poisson’s ratio, v 12 , and the laminate coefficient of mutual influence, rj 12 i 


respectively 
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Now, if the laminate strain e?, is applied to the individual layers, each layer will deform in 


accordance with its characteristic elastic properties 
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The stresses 622 and O12 develop because the strains in equation ( 3 . 13 ) are required to match 


the strains in equation ( 3 . 1 1). Expressions for the stresses arc developed by writing the strain in 


each layer as a combination of the individual ply strains e^, and terms, where the 


terms are required to force the total strains in each ply to match the laminate strains 




(3.15a) 



22 


t° 72 = e^ + Se$ (3.15b) 

= + «Y® (3.15c) 

Substituting into equation (3.10) gives 

<Ti? = Efi ) e?i +Q 12 &© +Qi65yf^ (3.16a) 

O 22 = Q22^ e S +Q265Y?^ (3.16b) 

on = Q^Seg + (3. 16c) 

From equations (3.11) and (3.13), 

= ( E 22 - ) = e ?i ( — v i 2 +v ^ ) (3.17a) 

= ( 7?2 - Y g ) = e?. (ti, 2 .i-Tl^i ) (3.17b) 


Combining equations (3.16) and (3.17) then gives the laminae stresses in terms of the laminae 
stiffnesses, the applied axial strain, and the mismatch in Poisson’s ratio and coefficient of mutual 
influence of the laminae and the laminate: 

du = [Efti + Q12 (vf^-Vi 2 ) + Q16 Cni 2 ,i ~ “nf&i )1 E 11 

O22 = [Q22 ( v ?2? _ V12) + Q26 (fll2.1 _T l?^,l)] E ?l ( 3 . 18 ) 

0\2 = [Q26 ( V ?2? -v 12) + Q66( T ll2.1 T lf^,l ) 3 E ?1 

From these equations it is seen that when the material constants V 12 and rin.i of the individual 
plies are identical to those of the laminate, the in-plane stresses 022 and <5 1*2 are zero. On the 
other hand, when there is a difference in the material constants, the in-plane stresses will in 
general be non-zero. As shown subsequently in Section 3.3, the magnitudes and signs of the in- 
plane stresses, and the laminate stacking sequence, have a direct influence on the magnitudes of 
the interlaminar stress components. 

3.2.3 Classical Lamination Theory Summary 


The classical lamination theory solution is approximate and in general only satisfies the 
equilibrium equations and edge boundary conditions in a through-thickness average sense. On a 
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point by point basis, however, the differential equations of equilibrium are not satisfied by the 
classical lamination plate theory solution in edge regions where there is a transverse gradient in 
the in-plane stress field, because the out-of-plane stresses are assumed to be zero. Further the 
surface tractions arc not zero as required by the exact elasticity equations. Thus, the stress field 
needs to be refined in regions near the boundaries. The next section details the formulation of 
the refined solution developed in the present investigation. 

3.3 Boundary Layer Stress Solution 

3.3.1 Problem Formulation 

The solution is developed by recalling that the applied loading and hence the stresses and 
strains are independent of the the axial coordinate. The analysis may then be restricted to any 
y-z cross section. To take into account warping of the cross section, induced by the presence of 
off-axis layers, the generalized plane deformation assumption, with orthogonal displacement 
components u,v, and w is employed/ 60 ^ As shown in Figure 3.1, the free edge is defined by y = 
0. We assume that the classical lamination theory solution has been obtained and concentrate on 
calculating stresses in the boundary layer region. The classical lamination theory solution is 

00 (k) 

valid in the interior but predicts non-zero stresses a 12 and O 22 at the free edge. This defect in 
the satisfaction of the free edge boundary conditions is corrected by a refined approximate 
solution that assumes the total stresses in each layer to be a combination of the classical 
lamination theory stresses plus an additional contribution to the stress field which is negligible 
outside the boundary layer region. The in-plane components of the boundary layer terms 
evaluated at y=0 are taken as the negative of the classical lamination theory values so that the 
free edge conditions are satisfied. Thus we assume 
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o0°(y,z) = ^(y.z) + Sy ’(y.z) ( 3 - 19 ) 

where: 

a(j^ (y ,z) is the total stress in the kth ply 
00 

dy (y ,z) is the clt stress solution in the kth ply 
00 

6;j (y,z) is the local solution for the kth ply 


Because of the difficulties and inefficiencies encountered with solutions for the free edge 
stress field based on displacement formulations, a stress formulation is presented. Under the 
generalized plane deformation assumption, and in the absence of body forces, the equilibrium 
equations that must be satisfied in each layer have the form 
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The equilibrium equations (3.20) can be satisfied identically by expressing the stress components 
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in each layer in terms of the stress functions ,60] <j> and iy, such that 
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The sixth stress component On is determined from the compatibility equations (3.21a-c) and the 
strain-stress relations as subsequently described in Section 3.3.5. The stress functions and 
constitutive relations can then be used in the remaining compatibility equations to yield the 
following pair of coupled governing partial differential equations for the stress functions. 
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where, 1 * 2 , L 3 , and L 4 , are linear differential operators defined as: 
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Bj and B 2 are determined from the end conditions. 


(3.24) 


In addition to satisfying equations (3.23) for each layer, the solution must satisfy conditions 
on the external surfaces of the plate as well as satisfy the conditions of continuous tractions and 
displacements along interfaces between adjacent layers of the laminate. The traction free 
conditions at the edges y =0 imply 
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ofi = CT 12 +Oi 2 =0 

+ <*22 = 0 (for y = 0) (3.25) 

o$ = aa = 0 

The traction free conditions on the top and bottom surfaces require 

a, 3=0 

CT23=0 (for X 3 = ±H) (3.26) 

ct 33 =0 

For perfect bonding between layers, the continuity of displacements and tractions at the layer 
interfaces impose the six additional conditions on the solution in each layer 


°? 3 +1) (y.t (k+1) ) = oJV(y.O) (3.27a) 

<T? 3 +1) (y,t (k+ ») = CT^(y,0) (3.27b) 

CT? 3 +,> (y,t <k+1) ) = a$(y,0) (3.27c) 

u (k+ ,) (yt< k+1 )) = u®(y,0) (3.27d) 

v (k+l)^y((k+l)j = v W(y t 0) (3.27e) 

wtk+n^y^+O) _ w W(yO) (3.270 


where the superscripts k and k+1 designate the ply above and below the interface, respectively. 
We have also assumed that classical lamination theory stresses are recovered in the interior. 
This assumption implies the additional conditions 

limCT i j(y,z) <|l) = 0 (3.28) 

y-*“ 

Finally at the ends of the laminate, the kinematic conditions w = eoXj or w = (koX 3 )X| are 
imposed, where to and Kb arc applied axial strain and curvature, respectively. 

3.3.2 Stress Assumptions 

An approximate solution to equations (3.23) subject to the boundary conditions specified in 
equations (3.25-3.28) is obtained by choosing a stress field that exactly satisfies the differential 
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equilibrium equations (3.20a-3.20c), and the stress boundary conditions (3.25-3.28). For clarity, 
explicit expressions for stresses rather than stress functions are presented. Results from previous 
investigations are used to guide the selection of appropriate stress forms. The previous analyses 
show large through-thickness gradients as well as large in-plane gradients in the stress field in 
regions near the free edge and in interfacial regions. In the present investigation, these 
perturbations in the classical lamination theory stress predictions are approximated by assuming 
two physical effects contribute to the stress field. The first effect is that represented by the 
Kassapoglou and Lagace (KL) solution. In this solution, the in-plane stresses predicted by 
classical lamination theory are used to formulate expressions for the out-of-plane stresses. As 
previously shown in Section 3.1 the classical lamination theory stresses develop because of the 
mismatch in engineering properties between the laminate and the individual laminae. The stress 
field obtained by adding the KL refinement and the classical lamination theory solution exactly 
satisfies the differential equations of equilibrium, the stress free boundary conditions, and the 
traction continuity conditions. Thus this contribution is referred to as the global mismatch or 
global equilibrium effect. The solution obtained from global mismatch considerations, however, 
assumes the laminate behavior to be qualitatively the same throughout the thickness of the 
layers, and consequently does not capture the large through thickness gradients in the layer stress 
fields near the interfacial surfaces. Further, the KL method predicts zero stresses at some 
interfaces where other analysis methods predict stresses of considerable magnitude. The second 
effect included in the present formulation is intended to relax some of the constraints imposed by 
the KL technique and provide an improved approximate theory. This effect is more local in 
nature than the KL contribution to the stress field, and relates the stresses developing near 
interfacial planes to the relative displacements of adjacent layers, arising because of the 
mismatch in engineering properties between these layers. Two material property mismatches are 
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considered; a coefficient of mutual influence mismatch, and a Poisson’s ratio mismatch. 
Herakovich^ 8 * has shown that these properties are the most important to consider, with regard to 
interlaminar stress development, for the mechanical loading problem. Poisson’s ratio is defined 
as 


^22 



(3.29) 


and the coefficient of mutual influence is defined as 
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where Sjj are the previously defined compliance coefficients in laminate coordinates. The 
coefficient of mutual influence mismatch is primarily responsible for the development of and 
0 \ 2 t while the Poisson’s ratio mismatch is primarily responsible for development of O 22 . O 23 . 
and 033 . Since global equilibrium and the free edge boundary conditions are satisfied by the KL 
solution, self equilibrating forms are assumed for the local mismatch contributions to the stress 
field. 


The stress components for each layer are chosen as product functions of the thickness 
variable z, and the in-plane variable y. This leads to the following functional form for each of 
the effects in the k 1 * 1 ply: 

a) global mismatch 

»£,<**) -fjPWgjPfr) i= 1,2.3. j = 2,3 (3.31) 

b) mismatch in coefficient of mutual influence 

5 S n ,(y. z ) = h? ) (y)lS t) (z) i=l. j = 2,3 (3.32) 


c) mismatch in Poisson’s ratio 
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Oi^/y.z) = m^CyJn^z) + mf , (y)pj ) (z) i = 2,3, j = 2,3 (3.33) 


The stress field 6(y,z) in the kth ply is then the sum of (a), (b) and (c) 

(k) 

®U (y- 7 -) = f? > (y)gf ) (z) + h^'Cy)! J ) (z) + mW (y)nf’(z) + mj ) (y)p§ t> (z) i = 1,2,3, j = 2,3 (3.34) 

Substituting equation (3.31) into the differential equations of equilibrium yields the following 
system of ordinary differential equations: 
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Similar equations are obtained relating the functions in equations (3.32) and (3.34) If equations 
(3.32) and (3.33) are substituted into the differential equilibrium equations, relationships like 
equations (3.35) and (3.36) are obtained for the h, r s, and m^s, and expressions similar to 
equations (3.37) and (3.38) are obtained relating the lys, nyS, and PyS. Note that equation (3.35) 
uncouples from equation (3.36) and equation (3.37) uncouples from equation (3.38). This 
decoupling was indicated previously in equations (3.22). Thus only four functions have to be 
assumed for each contribution; two in-plane variations and two through-thickness variations. 
The remaining functions are determined from the conditions (3.35-3.38). 

3.3.3 Global Mismatch (KL Solution) 


Following the KL solution, the first refinement to the stress field is made by assuming 
approximate expressions for the out-of-plane stresses, based upon the in-planc stresses. The in- 
plane stresses are taken to have the form 
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«$*<y.*) = f$(y)gfi?(z) (3.39a) 

«£p,(y.*)“fg(y)g©<*> (3.39b) 

where f 5 ^(y) and ffe*(y) are unknown functions of y and g?^(z) and g^Cz) are unknown 
functions of z. The form of the functions g$ (z) and g$ (z) is determined by imposing the free 
edge conditions (3.25) at y = 0 

^(y = 0,z) = -dl k 2 ) (z) (3.40a) 

^(y = 0,z) = -a^(z) (3.40b) 

where 6n ( z ) and 6 12 GO are the in-plane stress components predicted by classical lamination 
theory. For a general symmetric or unsymmetric laminate subject to a uniform extension or 
bending load, classical lamination theory predicts stresses that vary at most linearly through the 
thickness of each ply. Therefore, we have 


g$(z) = B?>z + B?> 

(3.41) 

g^(z) = Bf ) z + B?> 

(3.42) 

-00 -00 
012,-012. 

1 t oo 

(3.43a) 

B ?> = a£ 

(3.43b) 

-00 -00 
3 t oo 

(3.43c) 

B ?> = &£, 

(3.43d) 


and the t and b subscripts denote the top and bottom of the kth ply, respectively, and t (k) is the 
thickness of the k '* 1 ply. These expressions are then substituted into the equilibrium equations 
( 3 . 37 ) and (3.38) and integrated with respect to z to obtain expressions for the through-thickness 
variations of the interlaminar shear components o^Cy.z) and 023 ^ (y.z). After imposing the 
stress free conditions at the top and bottom surfaces of the laminate, the interlaminar shears may 


be written as 
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o?sL(y.z) = fiV(y)[B? t) z 2 /2 + B? , z + B?>] (3.44) 

^(y.z) = f^(y)[Bf } z 2 /2 + B?>z + B?>] (3.45) 


In a similar manner, the interlaminar normal stress is obtained from equations (3.38) and 
(3.45) 


Onjy.z) = f^(y)[B? ) z 3 /6 + Bfz 2 /l + Bg°z + B?>] (3.46) 


Constants B^, B^, and ,B^ are determined from the interfacial traction continuity conditions. 
Starting at the bottom free surface and working up, they have the form 


B§°= £ j^BP(t®) 2 /2 + B^pt^j k= 1,N— 1 

Bg° =£ JbSP< 1®) 2 /2 + Blpt^j 


k= l.N-l 


N 

b?> = 2 

j=k+i 


B^(t^) 3 /6 + B^t®) 2 # + ^ ) (t (j) ) (2) /2 + B$H W ] £ t <m) 


(3.47) 

(3.48) 

k= l.N-l (3.49) 


and B^, B^ and B^ are all equal to zero. 

00 

The stress field o 1J(t) (y,z) is now expressed in terms of the 2N unknown functions f^V(y) and 
f^(y), where N is the number of layers in the laminate. If f$(y) and f^*(y) are expressed in 
terms of quantities at the interfaces, so that interfacial continuity of tractions is guaranteed, a 
general solution for these functions may be obtained by invoking the principle of minimum 
complementary energy. According to this principle, out of all possible stress fields o^y.z), 
that satisfy both equilibrium and the stress boundary conditions, the one that represents the 
actual equilibrium state, is the one that minimizes the laminate complementary energy. 171 1 The 
complementary energy is defined as the strain energy of the laminate minus the external work 
done on the portion of the laminate where the displacements are prescribed 


nc-in^^xJJJnansHa} 

k=l k=l v w 


l (k)dVk -JJ[ {T ) T {u ldA 


(3.50) 
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Using the assumed forms for the stresses in equation (3.50) an expression for the complementary 
energy in terms of the arbitrary functions fjj^y) is obtained. Since we have assumed stresses to 
be independent of the axial coordinate Xj, and have assumed explicit functions for the through- 
thickness variations of the stresses, the z and x integrations can be carried out, reducing the 
volume integral in equation (3.50) to a line integral in y. Taking the first variation of the 
simplified integral and equating it to zero yields a system of 2N, constant coefficient, ordinary 
differential equations for the functions f$(y) and f$(y). These equations and the homogeneous 
boundary conditions along the free edge define an eigenvalue problem whose solutions are 
exponential functions of y. The complete solution for any of the functions f^ and is then 
obtained as a combination of a particular solution and a linear combination of the eigenfunctions 
for the homogeneous solution. 

A general eigenfunction expansion solution of this type is favorable for a few reasons. First, 
since separate functions are assumed for each ply, the stresses in individual plies can decay at 
different rates. Second, as Pagano has shown, the accuracy of the stress field predictions will 
improve significantly when the number of sublayers used to model a layer increases, just as 
results improve with finite element models when the mesh is refined. For the same through- 
thickness discretization, a model of this type will provide more accurate results than finite 
element models, since continuous functions are used to describe the y variation in the stress field 
rather than discretized functions as are used in the finite element models. These advantages, 
however, come at the high cost of solution inefficiency, which is precisely what we are trying to 
avoid. Run times are a function of the number of layers in the laminate and can become 
prohibitive for laminates with a large number of plies. Also, as noted earlier, numerical 
constraints have been shown to limit the number of plies in a laminate that can analyzed by this 
approach.^ 51 * 
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Considering the efficiency problems associated with the generalized eigenfunction 
expansion technique outlined above, KL took a slightly different approach and made some 
assumptions at the outset on the form of the functions f^(y) and f?^(y). These assumptions 
simplified the energy expression and subsequently reduced the system of equations that had to be 
solved. First, based upon the fact that the integrals of the in-plane stresses a $ and through 
the laminate thickness are zero at any y location, they assumed that ff^Cy) and f©(y) be the 
same in all layers 

^’(y.z) = a(z) k fi 2 (y), i=l,2 (3.51) 

By taking in-plane stresses of this form, the analysis greatly simplifies since the number of 
unknown functions is reduced from 2N to two. The limitation imposed by this technique, 
however, is that the stresses in all plies are forced to decay at the same rate. KL also went one 
step further, and assumed explicit sums of exponentials, expressed in terms of unknown decay 
parameters, for the in-plane variations f1z*(y) and f^(y). The interlaminar functions were 
derived from the in-plane functions and the differential equations of equilibrium. Minimization 
of the laminate complementary energy then resulted in a system of two non-linear algebraic 
equations for determining the unknown parameters. 

KL used exponential functions to describe the in-plane variation in the stresses because they 
provide for the necessarily rapid decay of the boundary layer stress components with distance 
from the free edge. Explicit sums of exponential functions are used in this study as well because 
the KL predictions generally show good trends, and as mentioned above, the Euler Lagrange 
equations obtained from the more general approach arc solved by a series of exponential terms. 

KL determined the necessary forms of the in-plane functions by enforcing the requirements 
of the Force-Balance Method.^ The Force-Balance method will also be used in deducing the 
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required forms of the in-plane variations associated with local mismatch effects. The Force- 
Balance Method is basically a statement of overall force and moment equilibrium, applied to a 
section of the laminate large enough to satisfy the assumptions of material homogeneity. If a 
rectangular volume element is taken with its X 2 faces at the laminate center plane and stress free 
edge, and with the X 3 + face corresponding to the top free surface, as shown in Figure 3.1, six 
equations are derived from the force and moment equilibrium conditions.' 20 ' 

Force Equilibrium: 


£Fx.: 

Ja 13 dy + J a 12 dz = 0 
3 - r 

(3.52a) 

11%: 

J + J = 0 

r r 

(3.52b) 

I*%: 

J a 33 dy = 0 

(3.52c) 


3 


Moment Equilibrium: 

£M Xl : J cr 22 zdz + J a 33 ydy = 0 (3.52d) 

r 2 r 

£M Xl : Joi 3 dydz+ Ja 12 zdz = 0 (3.52e) 

1 * r 

E M x, : J a 12 dydz + J a n (b-y)dy = 0 (3.520 

r z 

The subscript on the integral indicates the face over which the integrations are taken. 

Equation (3.52c) implies that o$(y,z) is a couple. Therefore, o^Cy.z) must cross the y axis 
at least once, and at least two exponential terms arc required to represent this component of 
stress. KL assumed that one crossing of the y axis represents the lowest energy solution so they 
used a sum of two exponential functions to represent the transverse variations f$(y), f$(y). and 
f^(y). Similar considerations were used to determine f^(y) and f^(y). 
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Considering the above requirements on the form of fjj(y) and the free edge boundary 
conditions, KL obtained the following final expressions for the stresses in the kth ply: 


-00 

CT 22 fli ) - 

1 1 

1 

1 “ 

* N " 

"*, C 

k J J 

[b?>z + B?>] 

-00 i2 *0 

a33 (E) -4>2 Xi _i 

X 1 e-^ , -e‘ +iy ][B?>z 3 /64 

B?>z 2 /2 + Bg°z + B? 0 


02I =* 2 ^[ - c- X ‘^ y ] [ B?>z 2 /2 + B? >z + B?> 

aSL=«,c^[BPz 2 /2 + BPz + Bp] 


(3.53a) 

(3.53b) 

(3.53c) 

(3.53d) 

(3.53c) 


where <|>i, <t» 2 , and Xj are unknown decay constants determined by minimizing the 
complementary energy. KL assumed 4>i = <j> 2 so only two unknown parameters 4>i and X| are 
used in their solution. 


These stress distributions superposed on the classical lamination theory solution satisfy 
pointwise and global equilibrium, stress free boundary conditions, and interfacial traction 
continuity conditions. Also equations (3.53) show that for large y all of the a,jS are zero, so 
outside the boundary layer region the classical lamination theory solution is recovered. Further, 
these equations show that the interlaminar stresses predicted by the KL solution at an interface 
are proportional to B§^, and B^. The B’s contain stacking sequence information and are 
basically force ( B^\ B^) and moment (B^) resultants of the classical lamination theory 
stresses obtained by carrying out the through-thickness integrations in equations (3.52). 
Therefore, if the resultant force at an interface, or at any through-thickness location, z, is zero, 
the KL solution predicts identically zero interlaminar stresses at all points along the intcrfacial 
plane having this z-location. This form of assumption, however, is too restrictive. The overall 
equilibrium equations (3.52) only require that the integral of the interlaminar stresses be zero at 
any z-location where the through-thickness integration is zero, and not that the stresses 
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themselves be identically zero at all points along the interfacial plane. In fact, as previously 
mentioned, results from earlier investigations indicate large interlaminar stresses in interfacial 
regions where local mismatch of properties is present but interlaminar stresses are not required 
to satisfy global equilibrium. The next sections outline extension of the KL solution to include 
additional terms associated with this local property mismatch in the assumed stress field. 

3.3.4 Local Material Property Mismatch 

Stress assumptions to include local mismatch considerations are developed by defining the 
interlaminar shear stresses and the interlaminar normal stress at each interface and deriving the 
remaining stress components from these definitions using the differential equations of 
equilibrium. This approach has the advantage that the traction continuity conditions are satisfied 
by the form of the assumptions. The most general form of a definition of this type would have 
different functions at each interface, as mentioned in the previous section, allowing the stresses 
in individual plies to decay at different rates. This is a considerable relaxation of the constraint 
imposed by the KL solution that stresses in all plies decay at the same rate and should lead to 
increased accuracy in the stress predictions. The drawback of course, is a large reduction in 
solution efficiency since the number of unknowns is dependent on the number of plies in the 
laminate. The present formulation is a compromise between the general formulation and that of 
KL. The same functions are used at all interfaces to define the local mismatch contribution to 
the stress field. As in the KL solution, using the same function for all interfaces requires the 
decay of the local mismatch contributions to the stress field to be the same in all plies. However, 
the decay of the total stress, i.e., equilibrium contribution plus mismatch contribution, may differ 
in individual plies, if there is a local mismatch contribution in one ply but not in another. 


The interlaminar shear stresses arising at an interface are assumed to be proportional to the 
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mismatch in coefficient of mutual influence or mismatch in Poisson’s ratio between the ply 

above and below the interface of interest. The interlaminar normal stress is assumed to be 

proportional to the mismatch in Poisson’s ratio. The mismatch in coefficient of mutual influence 

is assumed to affect only the and 0 $ components of stress while the Poisson’s ratio 

mismatch is assumed to affect only the o^ 1 , o^j*, and components of stress. This 

assumption is exact for cross-ply laminates, where only and are present, but is 

approximate for more general laminates. In angle-ply laminates, for example, the classical 

00 

lamination theory stresses o 2 2 are all zero and there is no Poisson’s ratio mismatch between 
adjacent plies, but interlaminar stress components and a$, and transverse in-plane stress 
have been shown to develop. I 9 - 26 - 35 - 59 ! These stresses, therefore must result from the 
coupling between the o^, and components of stress and the o^?, o§^, and 0$ stress 
components. This coupling is apparent from the compatibility equations ( 3 . 23 ). For cross-ply 
laminates, the compatibility equations uncouple and and 0$ are identically zero. The 
previous predictions, however, have shown o o^, and to be an order of magnitude 
smaller than 0 $ and in angle-ply laminates so neglecting this coupling is a felt to be a 
reasonable assumption. 

The in-plane functions, h^y), m^Cy), are, as in the KL solution, assumed to be an explicit 
combination of exponential functions which are chosen such that the interlaminar stress 
components arising from local mismatch integrate to zero over y. For the m^(y) functions, the 
additional restriction that m^ (y) has a zero moment about the longitudinal axis is imposed. 
Thus, the mismatch effect permits non-zero stress contributions, but these stress contributions do 
not alter the global force and moment equilibrium established through the fjj^y) and gjj^(z) 
functions of the KL solution. The through thickness variations lfj^z), n^(z), and p{j^(z) are 
polynomial functions chosen so that the stresses resulting from mismatch decay with distance 
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from an interface. In order to keep the number of unknowns in the solution to a manageable 
number, the through-thickness decay lengths of the mismatch effects are established a priori and 
set equal to the thickness I*** of the individual layers in the laminate. A more general 
formulation would allow the decay lengths to be variable and left as unknowns to be determined 
by minimizing the laminate complementary energy, but once again at the expense of increased 
computational cost. The effect of variable through-thickness decay lengths, obtained by 
dividing a layer into sublayers, is briefly addressed in Section 4.1.1. 

Mismatch in Coefficient of Mutual Influence. The interlaminar shear stress 0^3* at an 
interface, associated with the mismatch in coefficient of mutual influence, is assumed to be 
proportional to the mismatch in r| 12 .i of the two plies adjacent to that interface. The 
proportionality constant is a product of an unknown constant Aj, determined by minimizing the 
laminate complementary energy, and the applied axial strain ejj. The same constant is used for 
all interfaces. 

The mismatch in coefficient of mutual influence at the two interfaces bounding the k* layer 
are defined 

$■112,1 0 s * 0 = fliufk-l) - fli2,i(k) (3.54) 

8fli2.i(k»2) = Tli2.l(k) - rit2.i(k+l) 

with 5rii2,i(l.l) = 8r|i2,i(N,2) = 0. The layers above and below the kth layer are designated k-1 
and k+1, respectively, as indicated below. 

k-1 

k 

k+1 


The stresses in a generic layer are influenced by mismatch effects from the two adjoining 
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interfaces. Two quadratic functions are used to describe the through-thickness variation in shear 
stress within a layer. One function has the value of the mismatch at the top interface of the kth 
layer, and a value of zero at the bottom interface, while the other function has the value of the 
mismatch at the bottom interface of the kth layer, and a value of zero at the top interface. 
Therefore, the through-thickness function 1^ is assumed to have the form 


!$(z)= 8ni2.i(k,l)e t <k) z 2 /(t <k) ) 2 + 5Ti I2il (k,2)ef ) (l-z/t 0t) ) 2 A, 


(3.55) 


Equilibrium then requires that the function 1^ have the form 

• 5n 12 .i(k,l)e l w +5ri 1 2. 1 (k,2)eg t) j - 6Ti 12>1 (k,2)eg t) | A, (3.56) 




dz 


2z 

tOOi2 


where Et 00 and e§° are the strains at the top and bottom of the kth layer, respectively, and t 00 is 
the thickness of the kth layer. 


The associated in-plane variation h 13 (y) is chosen to be the same for all interfaces with the 
self-equilibrating shape shown in Figure 3.2. As previously mentioned, the same function is 
used for all interfaces to keep the number of unknowns to a minimum. To obtain the shape 
shown in Figure 3.2, a combination of two exponential functions is assumed 

hi 3 (y) = D 1 e'^ y + D 2 e- w (3.57) 


Integrating with respect to y and setting the result to zero, provides a relationship between D, 
and D 2 . 


Pi + P2 _ _ 

4*3 

Setting D ( = -1, hf^(y) has the self-equilibrating form 

h,3(y) = X 2 e-^ y -e- 4 ’ y 

It follows from equilibrium that h )2 (y) must have the form 


(3.58) 


(3.59) 
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Figure 3.2. Required form for h 13 (y) 

h, 2 <y) = J h I3 (y) dy = ± [-***♦" + e*"] + D 6 (3.60) 

where the constant of integration, D 6 , must be zero to satisfy the condition that the classical 
lamination theory solution is recovered in the interior. 

Mismatch in Poisson’s Ratio. Assumptions for stresses arising from a mismatch in 
Poisson’s ratio are derived in a similar manner to those developed for the coefficient of mutual 
influence mismatch. The mismatch in Poisson’s ratio at an interface results directly in an 
interlaminar shear stress, a^, and an interlaminar normal stress, Thus, there are two 
contributions to the assumed stress state arising from mismatch in Poisson’s ratio considerations. 
The first contribution is in the form of a direct assumption on the interlaminar shear stress at an 
interface, as was done with the coefficient of mutual influence mismatch. Interlaminar normal 
stress, and transverse in-plane stress, o$, are then derived from the differential equations of 
equilibrium. The second contribution is in the form of an assumption on o§9, with and 



41 


derived from equilibrium considerations. For both contributions, the stress assumed at an 
interface, either or a^, is assumed to be proportional to the mismatch in Poisson’s ratio, 
v 12 , of the plies adjacent to that interface. As with the coefficient of mutual influence mismatch, 
the proportionality constants are unknown and are determined by minimizing the laminate 
complementary energy. The same constants are used for all interfaces. 

The definition of the mismatch in Poisson’s ratio is similar to the definition for the 
coefficient of mutual influence mismatch 

Svi 2 (k>l) = Vi 2 (k— 1) — v 12 (k) (3.61a) 

5v 12 (k,2) = v 12 (k) - v 12 (k+l) (3.61b) 

with 5v 12 (l,l) = 5v 12 (N,2) = 0. 

Expressions for the through-thickness variations in the stresses in the k 1 ** ply are developed 
by considering a two ply, unsymmetric laminate. Figure 3.3 shows a section of the laminate 
with the shear stress a 23 (y,z) acting over a face parallel to the free edge and the interlaminar 
normal stress, a 33 (y,z), acting on the interfacial plane between layer (1) and layer (2). For layer 
equilibrium the integral 

«10 

J CT 23 (y,z)dz (3.62) 

o 

evaluated at y = y 0 must equal the integral 

yo 

Ja 33 (y,z)dy (3 63) 

o 

over the interface between the two layers. For laminate equilibrium the integral of o 23 (y,z) 
through the laminate thickness must be zero. (The latter restriction was not imposed when 
developing the through thickness assumptions for Oi 3 (y,z) in the previous section because equal 
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Figure 3.3. Laminate and Sublaminate Equilibrium 

and opposite shear stresses act on the ends of the laminate). Possible through-thickness 
functions for the stress component 623 (y,z) that satisfy these equilibrium requirements are 
shown in Figure 3.4. 

We begin by formulating the through-thickness functions n^(z) associated with the 
assumption on the interlaminar shear stress at an interface. In general there will be a non-zero 
strain e n at the (l)-( 2 ) interface, and consequently a non-zero shear stress C 23 there, since we 
are assuming o 2 3 (y,z) to be proportional to the axial strain and mismatch in Poisson’s ratio at 
the interface. In order to have a non-zero shear stress at the interface and still satisfy the 
requirement of overall laminate equilibrium, the through thickness variation of ^(y.z) must 
have the form shown in Case A of Figure 3.4. That is, n^(z) must integrate to zero over each 
layer thickness. If 623 (y.z) integrates to zero through the thickness of a layer at any location y, 
then c 33 (y,z) will be zero at all points y along all interfaces. 
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Figure 3.4. Through-thickness distributions of 023 (y,z) 

Two quadratic functions are employed to describe the through-thickness variation in shear 
stress within a layer as was used for the variation in the shear stress af 3 \y,z) in the previous 
section. An expression for the function associated with the mismatch at the bottom of a ply can 
be obtained by considering the first ply in Figure 3.4. We assume 

n & 2) (z) = Az 2 + Bz + C (3.64) 

Imposing the condition that n& 2) is proportional to the mismatch 5v 12 (k,2) at the bottom 
interface of the k 1 * 1 layer, and has a value of zero at the top interface gives 

CT 23 (y,z (1) = 0) = 6v 12 (l,2)e{, 1) (3.65a) 

CT 2 3(y.z (l) = t (1) ) = 0 (3. 6 5 b ) 

, . . . „ (k) 

and imposing the requirement that a 33 is zero at all interfaces as outlined above 
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a 33 (y,z (1) =t (1 >) = 0 (3.65c) 

a 33 (y,z (1 >=0) = 0 (3.65d) 


yields 

nS> k 3 ' 2) = 8v 12 (k,2)e£° ^ 3z 2 /(t°°) 2 - 4z/t (k) + 1 j A 2 (3.66) 

Similarly, considering ply 2 provides an expression for the function associated with the 
mismatch at the top of the k* ply 

n&° = 8v I2 (k, l)e,°° [ 3z 2 /(t w ) 2 - 2z/t (k) ] A 2 (3.67) 

The through-thickness variation n$(z) in each layer is then the sum of equations (3.66) and 
(3.67). 

n$(z)= [5v 12 (k,l)e l (k ^3z 2 /(t (k) ) 2 -2z/t (l) ] + 5v 12 (k,2)£g° [ 3z 2 (t w ) 2 - 4z/t w + lj] A 2 (3.68) 

where A 2 is an unknown constant to be determined by minimizing the laminate complementary 
energy. Equilibrium then requires that the functions n$-?(z) and n^(z) have the form 

n^(z) = [8v 12 (k,l)e«[6z/(t«) 2 -2/t^] +5v 12 (k,2)eg t >[6z/(t w ) 2 -4/t<“>]] A 2 (3.69) 


n$(z) = Sv 12 (k, l)e t w [ z 3 /(t«) 2 - z 2 /t°°] + 8v 12 (k,2)eg l) [ z 3 /(t w ) 2 - 2z 2 /t°° 


+ z 


(3.70) 


The corresponding in-plane function m 23 (y) is chosen to be self-equilibrating with the shape 
shown in Figure 3.5 To obtain a shape of this form, a combination of three exponential 
functions is required 

m23 (y) = D] eT* ,y + D 2 e^ y + D 3 c^ y (3.71) 


It follows from equilibrium that 


m22(y) = Jm 23 (y)dy = -^e^ ,y -— ^e~* <x,y _ ^^-e '* 4X<y+D 4 


(3.72) 
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Figure 3.5. Required form for m^Cy) 

m 33 (y) = dm ' 3(y) = -<t> 4 D,e^ y - <t> 4 X 3 D 2 e^ y - foX 4 D 3 e~**« y (3.73) 

dy 

where the integration constant D 4 must be zero so that m 2 2 (y) is zero for large y. We can now 
set Di to one, without loss of generality, and solve for D 2 and D 3 by imposing the free edge 
conditions. The free edge condition on m 23 (y) gives 

m 23 (y=0) — 1 + D 2 + D 3 = 0 (3.74) 


and the free edge condition on m 22 gives 

D 2 d 3 

m 22 (y=0) = -1/(J) 4 - — — - — r- 

<p 4 A. 3 <?4^4 

Dj and D 3 are then found to be 

_ Ml - K) X4(X 3 - 1) 

(X 4 - Xj) 3 (X 4 - X 3 ) 


(3.75) 


(3.76) 


The through-width functions associated with the Poisson’s ratio mismatch and the direct 
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assumption on the interlaminar shear stress then have the form 



»-4uy 


+ 


iiz. w ^ 

(X 4 - X3) 


+ 


(X 3 - 1 ) w 
(^4 ~ ^3) 


(3.77) 


m 23 


-♦«y , Ml ~ ^ r ^y , ~ 11 -~»4^y 

0 ^~h) (X4-X 3 ) 


m 33 - —<{>4 


eT* ,y + 


Ml -»«^r + 

(X4 - X3) 


- 1) . Ky 
(X4 — X. 3 ) 


(3.78) 


(3.79) 


Note these forms satisfy the self-equilibrating requirements outlined in the introduction to this 
section, provided classical lamination theory stresses are recovered in the interior region, i.e. 
provided the laminate is wide enough for the exponential terms to be zero in regions removed 
from the edge. Under these conditions, the self equilibrating requirements are identically met as 
a result of satisfying the boundary conditions. 

y=b 

J ni23(y)dy = m 22 (b) — 1 TI 22 (0) = 0 (3.80) 

y=0 

y 5 * 

J m 33 (y)dy = m 23 (b) - m^fO) = 0 (3.81) 

y=o 

y=b 

J m 33 (y)dy = m 22 (b) - m 22 ( 0 ) = 0 (3.82) 

y=0 


To develop expressions for the direct assumption on interlaminar normal stress we again 
consider a two layer unsymmetric laminate. We assume 033(y,z) varies cubically through a 
layer thickness so the interlaminar shear stress 023(y,z) will again vary quadratically. Recall 
that for laminate equilibrium and layer equilibrium to be satisfied, a23(y,z) must be distributed 
through the laminate thickness with a form like that shown in either Figure 3 . 4 a or 3 . 4 b. Since it 
is desired to have a non-zero interlaminar normal stress at the interfaces and since we have 
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already explicitly assumed 023 (y,z) at the interfaces, a^Cy.z) associated with the normal stress 
assumption will have the form shown in Figure 3.4b. Using this requirement, and proceeding as 
was done for the n^ function, the following through thickness variation p^V for the interlaminar 
normal stress in the k 1 * 1 layer is obtained 

p$(z) = [sv^fk.Oe^f— 2z 3 /(t ck) ) 3 + 3z 2 /(t (k) )j +8v 12 o (k,2)eg t) [2z 3 /(t <k) ) 3 - BzVft ^) 2 + 1 A 3 (3.83) 
Equilibrium then requires the functions pi$(z) and p^(z) have the form 

pgftz) = [ 6v 12n (k, l)e,° t> — 6z 2 /(t (k) ) 3 + 6z/(t°°) 2 ] + 8 v 12b (k,2)e£° J 6z 2 /(t (k) ) 3 - 6z/(t (i) ) 2 + 1 A 3 (3.84) 

p^(z) = [ 5v 12o (k, ljef) [-12z/(t«) 3 + 6/(t (i) ) 2 j + 5v I2 .(k,2)eg l) [ 12z/(t«) 3 - 6/(t w ) 2 ]] A 3 (3.85) 

The in equations (3.83-3.85) are the equal to the 8 vj 2 in equations (3.68-3.70) for layers 
above the midplane and equal to their negative for layers below the midplane, since the 
interlaminar normal stresses are symmetric with respect to the midplane for the case of uniform 
axial extension of a symmetric laminate, and are antisymmetric for the case of uniform bending. 
The through-width assumptions associated with equations (3.83-3.85) are of the same form as 
equations (3.77-3.79) but divided by <|> to give dimensionally consistent stress expressions. 

3.3.5 Total Stress Assumptions 

The final form of the stress assumptions in each layer of the laminate can now be obtained by 
combining the classical lamination theory solution and the stress field refinements obtained from 
global and local material mismatch considerations in accordance with equation (3.34). The final 
forms of the stress assumptions including all contributions arc then 
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The assumption that stresses are not dependent on the longitudinal coordinate X[ caused aft* 
to drop out of the equilibrium equations. The strain compatibility equations and strain-stress 
equations are used to determine aft*, as suggested in Reference [68]. It follows from equations 
(3.21a-c) that eft* is a linear function of y and z: 

eft) = A^y + B^z + C® (3.91) 

Now eft* can be expressed in terms of the stress field using the constitutive equations (3.1): 

eft) = Sn oft) + s£o» + S$og + = A w y + B«z + C® (3.92) 

For large y, the left hand side of equation (3.92) is independent of y because the classical 
lamination theory solution is recovered away from the free edge. This implies A^* = 0. Solving 
for aft* then gives: 

aft) = Kft)z + K?) - + S^ofgl (3.93) 

Sn L J 

The constants Kft* and Kft* are obtained by matching the solution at large y with the solution 
given by classical lamination theory. 

limaft) = Bft) Z + Bft/ = Kft)?. + Kft> - -=Jg-|s?2( B ? >z + B S° ) + Bft’z + Bft’ )] (3.94) 


We therefore obtain 



50 



r =« 

^OO 1 


r ^oc) 

=00 1 

aft> = 

BP + Ig-BP + ^BV 

Z + 

+ =£-B?> + 


^11 

^11 


^11 

11 


=ir s n°S + S “oH + ► 

S 11 l 


or more simply 


= fiCT 
^11 


i(k)-W — Ck) _(k) -<k)_Q0 

^llGn +J> 12 a22 +^16^12 


Si 2 cr^ + Si 3 a^ + S 16 a^ , 


(3.96a) 


where 


On= B^z + Bft? 
a^ = Bf>z + B?) 
af^Bfz + B?’ 


(3.96b) 

(3.96c) 

(3.96d) 


The stresses are now expressed in terms of the unknown decay parameters (fo and Xj and the 
proportionality constants A k , through equations (3.86-3.90) and equation 3.96. The condition of 
minimum complementary energy is used to determine these constants. 


3.3.6 Complementary Energy Minimization 

The complementary energy can be expressed as a summation over the individual plies as 

n c = £n?>= £}J J [ [a) T [S](aJ ] (k >dV k -J|{T) T {G)dA ( 3 . 97 ) 

k=l k=l V w S. 

where V k is the volume of the k* ply, S u is the portion of the boundary over which 
displacements are prescribed, { u ) are the prescribed displacements, {T) t are the associated 
tractions on the displacement boundary S u and [S]' k) is the compliance matrix of the k* layer 
with respect to the laminate coordinate system. 

In order to evaluate the energy expression in equation (3.97) the prescribed displacement { u } 
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must be determined. Since the displacement un in the laminate docs not vary with y over the 
ends where the displacements are prescribed, Un may be obtained by integrating the strain e u at 
the center of the laminate. That is 

u n = J e ndX] = J[SntTn + S 12 CT 22 + Si60i2]dX] (3.98) 

By symmetry, un =0, at X] = 0, so at X] = L, the prescribed displacement u lt is given by 

u n = (Si 1 Oi x +$ 12^22 +Si6^i2)L (3.99) 


The volume integration with respect to X] will yield a factor of 2L multiplying that term, since 
the stresses and compliances are not dependent on the longitudinal coordinate Xi and, therefore, 
it is only necessary to compute the complementary energy per unit length of the laminate. 
Additionally, for thin laminates 1721 only half of the laminate, (0<X 2 <b), needs to be considered. 
Further, since explicit piecewise continuous functions have been assumed for the through- 
thickness variations, the z integration is performed by summing the individual integrals in z^\ 
Making these simplifications and substituting for [S] , {u], and aft* we obtain the following 
expanded expression for the complementary energy in the laminate 


n c = Ejj f 

k-l V»L 


<*/ 


+ R 


33' 


a?/ a?/ ct^/ of/ 

'+R 6 6-f- + R44^-+R55^f- 


(3.100) 


+ + R260^CT?^ + R 3 60$a$ 


■= Si2<t 5$ + S 16 a^j ^S n CT n + S[2<T22 + S 16 a I2 j 


dV, 


00 


where the Ry are as previously defined in equation (3.24), and terms that arc not functions of the 
unknown parameters have been omitted. Once the integrations have been performed, the 
complementary energy n c may be written as 
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n c =C 1 (4. i> ^) + A 1 C 2 (4» i ,X,) + A 2 C 3 ((t> i ,^) + A3C4(«l>i,X i ) + A 1 A 2 C5(«t) 1 .X i ) (3.101) 

+ A 2 A 3 C 6 (<)) i ,X i ) + A] A 3 C7(4»j,^) + A]Cj(4>i,Xi) + A 2 C<((<|>j,X;) + A 3 Cio(4 ) i>^i) 

where the C m , (m = 1,10) are polynomial functions of 4>i and X, multiplied by constants d k , and 
are expanded in Appendix A. The polynomials represent the result of the integration in the 
transverse direction y, and the d k represent the integration through the thickness of the z 
variations of the stresses. Note, that the expressions for the transverse integrations presented in 
Appendix A, are not exact expressions for the definite integrals in equation (3.100), but are 
based on the assumption that 

e - * b = e~* |Xjb = 0 (3.102) 

This assumption implies that (fcb is very large, that is the laminate is wide relative to its 
thickness. This places a restriction on the geometry of laminates that can be analyzed using the 
present formulation, but the severity of this restriction cannot be ascertained until the parameters 
<j>i and Xj have been determined. 

The variation of n c with respect to the unknown parameters <|>j, Xj, and A k can now be taken 
and set equal to zero. This results in a system of fourteen coupled, non-linear algebraic 
equations for the unknowns written symbolically as 

^ = 0 (1=1,5) (3.103a) 

d<t>i 

-^- = 0 0=1.6) (3.103b) 

dAj 

an c 

— ^ = 0 (k=l,3) (3.103c) 

oAjj 

The form of these equations is shown in Appendix A. 

3.3.7 Solution of System of Equations 

The problem has now been reduced to solution of the systems of equations (3.103). Since 

these equations arc non-linear, root finding proceeds by iteration. Basically, two numerical 
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approaches can be employed for finding the roots. Since the equations represent the gradient of 
the complementary energy, one approach is to find the minimum of the energy function 
(objective function) using a general numerical unconstrained optimization technique. Another 
approach that can be taken is to solve the simultaneous system using Newton’s Method. 
Theoretically, the minimization technique should be more efficient, since the search for a 
minimum can basically be reduced to a one dimensional problem, that is, a minimum can be 
found by moving "downhill" on a single surface . 1731 There is not an analogous procedure for 
finding a multi-dimensional root. However, both methods were tried in the present 
investigation, and greater success was experienced using Newton’s method. 

Both Newton’s method and numerical optimization techniques are started by providing the 
iterative algorithm with an initial approximate guess to the solution. One difficulty with 
Newton’s method is that the solution may not converge from a given initial guess. Convergence 
is only guaranteed if the initial approximate solution is in the neighborhood of the solution. 
Another problem associated with the solution of non-linear equations is the possibility of 
multiple solutions, in our case corresponding to local minima of the energy function, requiring 
the iterative procedure to be initiated from several different starting values of the independent 
variables to ensure that the "best" solution is obtained. The "best" solution is the one 
corresponding to the lowest energy. In the present study, the initial approximate solutions were 
generated based upon the results of the KL solution. The KL solution can be solved very 
quickly and gives an indication of the magnitudes of the decay parameters <(>; and X,. Several 
starting points were generated by bracketing the KL solution and then incrementing the 
independent variables, <J>; and Xj within the "brackets”. Initial approximations for the 
proportionality constants A k were obtained by simultaneously solving the three equations 
(3.103c), at the starting values of <(>j and Lj. With the decay parameters known, the constants 
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Qn, (m=l, 10) can be evaluated, and equations (3.103c) are linear in the unknowns A k . Negative 
values of and Xj were discarded because they physically correspond to growing stresses with 
distance from the free edge, and complex roots were discarded because they in general lead to 
complex energy. 

3.3.8 Solution Implementation 

The equations presented in the previous section are incorporated in a FORTRAN program 
AAIS. The laminate configuration, i.e stacking sequence, and ply thicknesses, ply orientations, 
ply material properties, and loading conditions are required input. With this information, the 
classical lamination theory solution is obtained. 

Once the classical lamination theory solution has been obtained, the energy expressions may 
be formulated and solved for the unknown constants. As mentioned in the previous section, two 
approaches were used to obtain the constants. IMSL [741 routines were called to implement both 
approaches. For the optimization, IMSL routine DBCODH was used. DBCODH minimizes a 
function of N variables with simple bounds using a modified Newton method and a finite 
difference approximation to the Hessian. To solve equations (3.103) as a non-linear system, 
IMSL routine DNEQNF was used. DNEQNF uses a variation of Newton’s method and the 
finite-difference method to estimate the Jacobian. There are also IMSL routines available that 
require the user to provide an exact Jacobian. These methods were not used because of the 
length of calculations necessary to obtain the Jacobian matrix. 

The program output includes the classical lamination theory stresses, values for the unknown 
constants, the laminate strain energy, complementary energy and compliance constants. This 
information is then input to a postprocessor that can be used to calculate through thickness and 
interlaminar stress distributions at any desired location. Contributions from all effects are 
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delineated in the output files. 
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CHAPTER 4 


COMPARISON WITH PREVIOUS RESULTS 


In this chapter and in Chapter 5, the approximate theory presented in the previous sections is 
applied to the stress analysis of several finite width, straight free edge laminates subjected to 
uniform extension and bending. Because of the general nature of the formulation and its 
simplicity in terms of the number of unknown parameters, the applications are straight forward. 
Recall from the previous discussion, that for general laminates, the stresses in individual layers 
are expressed in terms of fourteen unspecified parameters, <(>; (i=l,5), (j=l,6), and A* (k=l,3) 

that are determined by minimizing the laminate complementary energy. For cross-ply and 
angle-ply laminates, in which some of the stress components are zero, or are assumed to be zero, 
fewer parameters must be determined. For cross-ply laminates the number of parameters 
reduces to ten, since a 12 and a i3 are identically zero, and in symmetric angle-ply laminates 
subjected to uniform extension the number reduces to three. The reduced number of unknown 
parameters for the angle-ply laminates is a reflection of the assumption made in the present 
formulation that is the only non-zero interlaminar stress. In unsymmetric angle-ply 
laminates and angle-ply laminates subjected to bending load, the number reduces to six. For a 
particular type of laminate, however, the number of unknown parameters is independent of the 
number of layers in the laminate, their material properties, and orientations. Obtaining classical 
lamination theory stresses in the plate interior is the only laminate dependent calculation which 
must be performed prior to determination of the interlaminar stresses. Consequently, 
computation times are linearly proportional to the number of layers in the laminate, and the 
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analysis can be applied to laminates with a large number of layers. 

As indicated in the previous chapter several simplying assumptions were made in 
formulating the present approximate theory in order to obtain a methodology with the efficiency 
characteristics outlined above. These assumptions will, to some extent, result in a decrease in 
accuracy of the response predicted by the present theory as related to that given by more 
complicated models. Therefore, to gain some confidence in the model, this chapter presents 
comparisons between predictions of the present theory and well known solutions available in the 
literature. First comparisons are presented for laminates subjected to extension loading. These 
comparisons are followed by results for laminates in uniform bending. 

4.1 Laminates in Uniform Extension 

In this section we compare the response predicted by the present method of analysis with that 
given by other investigators^ 35,59,7 ^ for finite width symmetric laminates subjected to uniform 
axial strain, en = 0 . 1 %. Specifically, four laminates are examined; [ 90 50 A) 5 o]s. [O50/9O50],, 
[±455 o]j. and [4550/— 1 1550/O50/9O50],. In all laminates, the layers have equal thickness h = 0.25 
in., and the thickness to width ratio of the laminate is taken to be one to four. Thus for the 
cross-ply and angle-ply laminates we have H = 2 h and b = 8 h, where a quasi-isotropic laminate 
has H = 4 h, b = 16 h. The clastic properties of each graphite-epoxy laminae are taken to be equal 
to those given in the early studies of interlaminar stresses^ 26 ! and are provided in Table 4 . 1 . 
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TABLE 4.1. Lamina Material Properties for Typical Graphite-Epoxy 


Typical Graphite-Epoxy Material Properties 

E, 

(msi) 

e 2 

(msi) 

G]2 

(msi) 

G23 

(msi) 

v 12 

v 2 3 

20.0 

2.1 

0.85 

0.85 

0.21 

0.21 


4.1.1 Angle-Ply Laminates 

The [±4550 ] s laminate has been studied by numerous investigators. In this section 
predictions given by the present theory for various stress components are compared with 
solutions obtained by Wang and ChoiJ 59 ^ using an eigenfunction expansion solution, and by 
Wang and Crossman, [3S| using a finite element analysis based upon constant strain, triangular 
elements. In all analyses the quasi-three dimensional assumption is made, that is, stresses and 
strains are assumed independent of the axial coordinate Xj. Before presenting results the issue 
of variable through-thickness decay of the local mismatch effects discussed in Section 3.3.4 is 
briefly addressed. The terminology and notations introduced in the following discussion will be 
used throughout the remainder of this text 

Recall from Section 3.3.4, that in order to keep the number of unknowns in the present 
solution to a manageable number, the through-thickness decay lengths of the local mismatch 
effects were established a priori and set equal to the thicknesses h (k) and h (k+,) of the layers 
adjacent to an interface where there is a mismatch in material properties. This specification of 
decay length is evident upon examination of equations (3.86-3.90). The extent in the thickness 
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direction of the local mismatch contributions can, however, be varied by representing each 
physical layer in the body as an assemblage of sublayers. For [0i/0 2 ] s or [0i/02lt laminates in 
which there is only one interfacial plane where there is a discontinuity in material properties, 
two sublayers in each layer with thicknesses ti and t 2 , as shown in Figure 4.1, are all that is 
required to assess the influence of the through-thickness decay length of local mismatch effects 
on the stress field predictions. For more general laminate configurations, where the stresses in a 
generic layer are influenced by mismatch effects from two adjoining interfaces, at most three 
sublayers are required. The results presented throughout this thesis were generated by either 
modeling each layer as one unit, i.e., no sublayers, or by dividing each layer into two sublayers 
with thicknesses tj and t 2 , as shown in Figure 4.1. Predictions obtained by representing each 
physical layer in the body as one unit are denoted by N=l, and results obtained by representing 
each layer as two sublayers are denoted N=2. 



Figure 4.1. Through Thickness Discretization 

The influence of the through-thickness discretization on the stress field predictions and 
laminate complementary energy has been studied for the [±45 J, laminate. As shown in Figure 
4.1, the decay length of the local mismatch effects will be equal to tj when two sublayers arc 
used to model a layer. The "optimum" values for tj and t 2 can be determined, for a particular 
laminate, by allowing them to vary in several applications of the AAIS program, and then 
selecting the set that corresponds to the minimum complementary energy. For the [±45], 
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laminate considered here, the values ti =0.11 in., and t 2 = 0.14 in. were obtained. The 
magnitude of a 13 at the intersection of the ±45 interface and the free edge has also been 
determined. Values are provided in Table 4.2 which shows the increase in the maximum 
magnitude o 13 , with decreasing decay length, down to the length corresponding to the minimum 
complementary energy. 

The through-thickness decay length not only affects the stress magnitudes, but also the 
distribution of stress. In order to illustrate this effect, two curves are presented in the subsequent 
figures (Figures 4.2 and 4.4) for the present theory predictions. The curves N = 1, as previously 
mentioned, correspond to modeling each layer as one unit, and the curves N = 2 were generated 
using the "optimum" thicknesses given above. 

The width distribution of the interlaminar shear stress o 13 at the ±45 interface predicted by 
the present technique and the analytical solution in [59] are shown in Figure 4.2. As can be seen, 
the present solution predictions for both N = 1 and N = 2 agree well with the previous solution. 
There is a small difference in the predictions at the intersection of the interface and the free edge. 
The elasticity solution of Wang and Choi^ predicts a stress singularity at this point. Similar 
behavior is displayed by displacement based finite element formulations, which predict 
increasing stresses in elements adjacent to the singular point, as the size of the elements near this 
location decreases. The present solution does not include a singularity in the formulation, but as 
previously discussed, the magnitude of ct ) 3 at the intersection of the interface and the free edge 
is a function of the through-thickness decay length of the local mismatch effect. 

Interlaminar stress distributions along the 45/-4S interface for and a 33 arc provided in 
Figure 4.3. For an angle-ply laminate, Poisson’s ratio mismatch between adjacent plies is zero, 
and the in-plane stress o 22 is zero, so the present theory predicts identically zero stresses a 23 and 
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TABLE 4.2. Increase in c 13 with Decrease in Decay Length 


t 2 (in) 

-Oi 3 (ksi) 

0.25 

1.43 

0.20 

1.52 

0.15 

1.64 

0.14 

1.67 

0.13 

1.70 

0.12 

1.74 

0.11 

1.78 


a 33 (see equations 3.89,3.90). The most pertinent observation to be made from this figure is that 
although the other solutions predict non-zero stresses O 23 and c 33 they are small compared with 
o 13 . The normal stress is less than 20% of o )3 and the shear stress O 23 is less than 5%. Note 
also that the cr 23 distribution predicted by Wang and Crossman^ 35 ' does not satisfy the traction 
free edge condition. Wang and Choi [59 l on the other hand, predict o 23 to be zero at the free 
edge, but their distribution for O 23 does not satisfy the transverse integral force equilibrium 
equation (3.52b), since a 23 obviously does not integrate to zero over X 2 . Differences are also 
noted in the interlaminar normal stress predictions. Wang and Crossman’s finite element 
solution predicts tensile stress at the free edge, but Wang and Choi’s elasticity solution predicts 
compressive stresses. 


Comparisons of through-thickness distributions of the interlaminar shear stress O 13 arc also 
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Figure 4.2. Comparison With Previous Results [59] for O13 at 45/-45 Interface in 
[+455 o/- 455 0 ] s Laminate (e n =0.1%) 

made. In Figure 4.4, distributions at X 2 /b = 1, and X 2 /b = 0.89 calculated by the present 
approach and Wang and Choi’s eigenfunction expansion solution, arc shown. The present 
solution for N=2 agrees fairly well with the elasticity solution at X 2 /b = 1. The major 
discrepancy in the results occurs at the intersection of the interface and the free edge, where the 
present analysis predicts finite maximum stress and the eigenfunction expansion solution 
becomes unbounded as X 2 — > b, and X 3 — » h. Away from the free edge, at X 2 /b = 0.89, 
generally good agreement between the two solutions is observed throughout the laminate 
thickness for both N=1 and N=2 of the present theory. 







X 2 /b 


(a) ct 2 3 Comparison 



X 2 /b 


(b) C 33 Comparison 

Figure 4.3. Comparison With Previous Results [35], [59] for o 2 3 and a 33 at 45/-45 Interface in 
[+45 5 oM5 5 o], Laminate (£ n =0.1%) 




-o 13 /e x (msi) 


(a) X 2 /b = 1.0 



-Oi 3 /e x (msi) 


(b) X 2 / b = 0.89 

Figure 4.4. Comparison With Previous Results [59] for Through-Thickness Oi 3 in 
[-r-45so/ — 455 q] s Laminate (£n =0.1%) 
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4.1.2 Cross-Ply Laminates 

Interlaminar stress comparisons between results predicted by the present theory and finite 
element analysis of Herakovich et al., |75] for [Oso/Wsolj and [905o/05ol» laminates are presented 
in this section. For these laminates, the in-plane shear stresses 012 and interlaminar shear 
stresses 013 are zero. The present solution results were obtained using N=l. 

Figure 4.5 shows the width dependence of the interlaminar shear stress O23 and normal stress 
O33, at the 0/90 interface of the two laminates. The present predictions compare fairly well with 
the finite element results, again differing mainly right at the intersection of the interface and the 
free edge. Note, the present theory satisfies the traction-free boundary condition exactly, and 
classical lamination theory stresses are recovered in the interior. The finite element results 
satisfy the traction-free boundary condition only approximately. There are also differences in 
the predictions of o 33 at the free edge with the most noteworthy being that the finite element 
method predicts tensile a 33 in the [90 4 /04] $ laminate where the present solution predicts a small 
compressive stress. The gradient of the 033 distribution is very steep near the free edge and the 
two solutions predict similar slopes. 

4.1.3 Quasi-Isotropic Laminates 

Quasi-isotropic laminates have been studied by a number of investigators using the finite 
element method. Here we consider specifically the laminate stacking sequence 
[455 o/-455o/05o/905o] s . Illustrative results for this laminate predicted by the present theory arc 
compared with Wang and Crossman’s^ 35 ^ finite element predictions. Each layer was represented 
as one unit (N=l) in obtaining the present results. 


Distributions of interlaminar normal stress O33 along the laminate midplane and the 45/-45 



°33 fell 
(msi) 



(a) (t 33 Comparison 


O23 fen 
(msi) 



(b) 023 Comparison 


Figure 4.5. Comparison With Previous Results [75] for a 2 3 and O33 at 90/0 Interface 
[9O4/O4L and [O4/9O4 ] s Laminates (Bn =0.1%) 




Figure 4 . 6 . Comparison With Previous Results [ 35 ] for 033 at Various Interfaces in 
[4550M550/O50/9O50]. Laminate (£n = 0 . 1 %) 
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Figure 4.7. Comparison With Previous Results [35] for Oi 3 at 45/-4S Interface in 
[455o/— 45 5o/0 5 o/90 5 o ] s Laminate (e u =0.1%) 

interface are shown in Figure 4.6. Clearly, both methods agree quite well at the midplane. 
Similar results were obtained at the interfaces X 3 = h. and X 3 = 2h. The predictions of the two 
methods at the 45/-45 interface (X 3 = 3h), however, differ significantly. The present analysis 
only shows the normal stress crossing the X 2 axis once, while the analysis in [35] predicts two 
crossings, with the stress reversing sign and becoming compressive near the free edge. Possible 
causes for this difference arc discussed in the next chapter. 

Figures 4.7 and 4.8 show comparisons of the interlaminar shear stress predictions, at the 
45/-45 (o 13 ) and 0/90 (o 23 ) interfaces, respectively. Again, for both cases, the present theory 
agrees quite well with the finite element results. 
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Figure 4.8. Comparison With Previous Results [35] for O23 at 0/90 Interface in 
[455 o/— 455 o/ 05 o/ 905 o] s Laminate (En =0.1%) 


4.2 Laminates in Bending 

As previously discussed, interlaminar stress calculations for laminates in bending have 
received relatively little attention in the literature. However, bending loads are common in 
practical applications. Salamon^ has presented results based upon a finite difference solution 
similar to that presented by Pipes and PaganoJ 26 ^ for a [05 o/ 905 o] s laminate with the material 
properties given in Table 4.1 and subjected to end moments about the transverse (X2) axis. 
Chan and Ochoa have presented results for several [O2/02 J* laminates under the same 
moment loading. Comparison is made here with the solution in [64], 

Comparison of the interlaminar normal stress along the 0/90 interface is shown in Figure 4.9. 
The laminate is loaded such that the maximum bending strain En =-0.1% is developed at the 
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top surface. Stresses arc normalized, as in [64], by the elementary bending stress (On ) that 
would develop in a undirectional laminate (0=0) at the X 3 location corresponding to the 0/90 
interface (dn = 10 ksi). As the figure shows, the agreement between solutions is not as good as 
in the axial extension case, but the trends are similar. The major discrepancy occurs at the free 
edge, where the finite difference solution predicts tensile stress about 1-1/2 times larger than the 
present solution. 
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CHAPTER 5 

APPLICATIONS AND DISCUSSION 


The comparisons presented in the previous chapter showed that the predictions of the present 
solution are, generally, in good agreement with predictions obtained from a variety of other 
solution methods for extension and bending loads in a variety of laminates. In this chapter, the 
analysis is applied to additional laminates to demonstrate the utility of the method as a design 
tool and to provide a brief study of the influence of laminate configuration and loading on the 
interlaminar stress state. Also, as was discussed in Chapter 3, the total stresses in the laminate 
arise from two physical mechanisms - local and global mismatch in material properties. In the 
subsequent discussion the relative significance of these contributions to the stress field is 
examined. 

Results in the form of through-thickness and interfacial stress distributions are presented for 
symmetric and unsymmetric cross-ply laminates and for symmetric angle-ply and quasi- 
isotropic laminates subjected to bending and extension loads. Through-width distributions are 
presented only at ply interfaces. Distributions for other locations, however, can be obtained 
easily. Comparisons are presented for results obtained using the present formulation, the KL 
solution, and finite element solutions, where available. As subsequently shown, the current 
solution predicts interlaminar stresses which are generally in close agreement with finite element 
results, and improves the KL solution primarily at interfaces where global equilibrium is 
satisfied by the lamination theory stresses, but local mismatch in material properties induces 
large interlaminar stresses. When studying the figures, recall that the magnitudes of stresses 
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predicted by the finite element models at the intersection of layer interfaces and the free edge 
(singular point) are a function of the fineness of the mesh near this point. Thus emphasis should 
be placed on comparison of entire stress distributions away from these points where finite 
element results are more accurate. 

Unless otherwise noted, the finite element results for extension of cross-ply laminates were 
generated using a program previously developed by Norwood and the finite element results 

f 'T'Tl 

for the angle-ply and quasi-isotropic laminates were generated using the program CLFE2D. 
Norwood’s program is based upon a full three-dimensional formulation, while CLFE2D makes 
the quasi-three dimensional assumption. Quarter symmetry models (i.e. one quarter of cross 
section modeled) were used in all analyses. 

For the symmetric laminates the axial strain loading was e, = 0.1%, and the axial curvature 
loading was K x = 0.1. For the unsymmetric laminates, it was more convenient to apply end 
loads rather than end strains and curvatures to obtain the classical lamination theory results. For 
these cases axial extension load N x and uniform bending load M x were applied. The material 
properties of a T300-5208 graphite epoxy used in the analyses are provided in Table 5.1. All 
plies were taken to have thickness t = 0.005 in. 

5.1 Cross-ply Laminates 

Cross-ply laminates are the simplest of the laminate configurations because there are no off 
axis plies. Consequently, the coefficient of mutual influence is zero in all plies and the 
interlaminar shear stress. On, is identically zero throughout the laminate. These laminates arc 
therefore studied to isolate the influence of Poisson’s ratio mismatch on the development of 
interlaminar shear stress o 2 3 and interlaminar normal stress a 33 . 
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TABLE 5.1. T300-5208 Graphite-Epoxy Material Properties 


T300-5208 Graphite-Epoxy Material Properties 

Ei 

(msi) 

e 2 

(msi) 

G )2 

(msi) 

G 23 

(msi) 

v 12 

v 23 

19.2 

1.56 

0.82 

0.52 

0.24 

0.49 


Results are presented for three laminates: [0 4 /90 4 ] s , [90 4 /0 4 ]5, and [0 4 /90 4 ] t , subjected to 
bending and extension loads. Classical lamination theory stresses for the symmetric and 
unsymmetric laminates for these two load cases are provided in Appendix C (Tables C.l and 
C.2). With the lamination theory stresses given, the unknown parameters in the assumed stress 
expressions have been determined by minimizing the laminate complementary energy. The 
values obtained for the present theory, along with the parameters <|>2 and X] required in the KL 
solution are given in Tables 5.2 (symmetric laminates) and 5.3 (unsymmetric laminates). The 
KL solution parameters are shown within parenthesis. Note, 4>i . fa, X -2 and do not appear in 
the tables since Oi 2 and O 13 are zero for these laminates. All results presented for the current 
method were obtained by representing each layer as one unit (N=l). The layers were all of equal 
thickness, h = 0.020 in. Analyses made with N=2, where each layer was divided into two 
sublayers of equal thickness, resulted in higher energies. The through-thickness decay length for 
N=1 then corresponds to the number of plies in the individual layers of the laminate. The decay 
lengths, denoted by h, arc also presented in Tables 5.2 and 5.3. 

The symmetric and unsymmetric laminates respond very differently to applied loads. In 
symmetric laminates the membrane-flexural coupling terms, B;j arc all zero, so that when these 




Solution Parameters For Symmetric Cross Ply Laminates 


Constant 

Extension 

Bending 

[O 4 / 9 O 4 ], 

[904/04]* 

[O 4 / 9 O 4 ], 

[ 9 O 4 /O 4 ]* 

<t>2 O/in) 


52.59 

(58.63) 


m 

<t>4 (1/in) 

28.30 


33.13 

45.15 

<t>5 (1/in) 

33.42 

36.39 

m 

45.73 

h 

I ' 

0.952 

(3.952) 

1.979 

(1.000) 

1.758 

(1.000) 

h 

3.817 

1.754 

8.025 

4.422 

X4 

11.49 

9.406 

m 

4.307 

*5 

2.175 

3.341 

3.326 

3.291 


1.905 

1.817 

1.543 

1.313 

a 2 xio -5 

(psi) 

-4.151 

-1.311 

-4.705 

-4.381 


A 3 XIO " 5 

-2.919 

3.988 

-6.104 

3.870 

(psi) 





h (in) 

.020 

.020 

.020 

.020 





















































Solution Parameters For Unsymmetric Cross Ply Laminate 


Constant 


Extension 


Bending 


<t >2 (1/in) 


68.72 

(201.4) 


70.33 

108.2 


<t) 4 (1/in) 


<t>5 (1/in) 




52.79 


89.25 


0.915 

( 1 . 000 ) 


54.73 


60.69 


2.660 

( 1 . 000 ) 




3.448 


2.831 


X 4 


3.632 


3.712 


*5 


1.312 


5.911 


^6 


A 2 X10" 5 

(psi) 


A 3 X10“ 5 

(psi) 


h (in) 


0.862 


-7.158 


5.074 


.020 


0.726 


-7.493 


-3.685 


.020 
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laminates are subjected to extensional loads they remain plane. The unsymmetric laminate, on 
the other hand, has non-zero Bn and B 22 and will deflect out of plane when extended. 
Comparison of stress distributions in symmetric and unsymmetric laminates then provides an 
indication of the effect of out-of-plane deflections on the interlaminar stresses. Norwood' 76 ' has 
shown that when the out of plane deflections are large - on the order of the laminate thickness - a 
nonlinear analysis which accounts for geometric coupling effects is required to accurately 
characterize the interlaminar stress response. In the present analysis, all out of plane deflections 
are assumed small so that a linear analysis is valid. 

5.1.1 Symmetric Laminates - Extensional Load 

Stress distributions for [0 4 /90 4 ], and [90 4 /0 4 ], laminates are provided in Figures 5. 1-5.8. 
Stresses determined by the finite element analysis, the KL solution, and the present solution, 
equations (3.88-3.90), are shown. Through-thickness distributions are provided for the top half 
of the laminate, with stresses plotted as a function of the normalized distance X 3 = X 3 /h from the 
laminate midplane, where h denotes one layer thickness. Interlaminar shear stress is anti- 
symmetric about the midplane and interlaminar normal stress is symmetric. Interfacial 
distributions are also shown, with stresses plotted as a function of the normalized distance X 2 /b 
from the laminate center, where b is the laminate half width. In all plots, stresses are normalized 
by the average far field stress N U /2H where Nn is the far field load obtained from lamination 
theory, and H is the laminate half thickness. The normalized stresses are denoted C33 and 623 . 

Figure 5.1 shows comparison of finite element results, the KL predictions and the present 
solution for the through-thickness distributions of the interlaminar normal stress O33 at 
X 2 = 0.999b for the two laminates. Evident from these figures is the significant improvement of 
the present theory over the KL solution, with the present solution predicting trends similar to 
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those displayed by the finite element results. In particular, note the asymmetry in response of 
the 0/90 and 90/0 laminates predicted by both the finite element analysis and the current model, 
while the KL solution, which is based entirely upon global equilibrium predicts close to 
symmetric response of these two laminates. That is the KL solution predicts the stress response 
in the two laminates to be basically mirror images of each other. Also significant from a design 
point of view, is failure of the KL solution to predict tensile interlaminar normal stress, at any 
location in the [90 4 /0 4 ] s laminate, since tensile interlaminar normal stresses are more 
detrimental to the integrity of laminated composite stracturcs than compressive interlaminar 
normal stresses. 

The difference in the stress response of the two laminates and the predictions of the three 
methods is further illustrated in Figure 5.2 which presents 033 distributions along the 0/90 
interface. Clearly, the present theory agrees quite well with the finite element results for both 
laminates, but the KL solution diverges near the free edge in the case of the [90 4 /0 4 ], laminate. 
Consistent results from all three solutions are obtained for this laminate only at sufficiently large 
distance from the free edge (e.g. X 2 = 0.95b). 

Similar comparative through-thickness and interfacial distributions for the interlaminar shear 
stress 023 are shown in Figures 5.3 and 5.4. Again, the results of the present solution compare 
more favorably with the finite element results than do the KL results. The most noticeable 
discrepancies arc observed in the through-thickness distributions. The relative difference in the 
maximum value of the shear stress obtained by finite elements and the present solution is 10%, 
for the 90/0 laminate, while the relative difference between finite elements and the KL solution 
is approximately 38%. Also note that the stress gradients in the thickness and width directions 
predicted by the present theory and finite clement method arc more severe than estimated by the 


KL solution. 
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Figure 5.1. Comparisons for o 33 at X 2 /b = 0.999 in [0 4 /90 4 ] s and [90 4 /0 4 ] s Laminates 
Extension 










X 3 /h 
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Figure 5.3. Comparisons for a 23 at X 2 /b = 0.993 in [O4/9O4], and [90 4 /0 4 ], Laminates 
Extension 
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The inability of the KL solution to predict the difference in behavior of the 0/90 and 90/0 
laminates described above suggests that this solution does not include enough degrees of 
freedom in the assumed stresses to accurately define the stress response, in some cases, and that 
"local effects" might be the cause of the asymmetry in stresses observed in these laminates. The 
relative influence of the local and global effects is illustrated in Figures 5.5 and 5.6, which 
delineate all of the contributions to the stress field incorporated in the present theory. The 
equilibrium contribution represents the first term in equations (3.89,3.90), the shear mismatch is 
the second term (multiplying A2) and the normal mismatch is the third term (multiplying A3). 
Note that although the KL solution and the equilibrium solution are identical in form, the KL 
solution curves in Figures 5.1 and 5.2 are not coincident with the equilibrium contribution plots 
in Figures 5.5 and 5.6 since different <(>'s and X's were employed to obtain the two distributions 
(see Table 5.2). 

Figure 5.5 provides through- thickness distributions for the interlaminar normal stress 033 
near the free edge. Equilibrium and mismatch effects contribute throughout both laminate 
thicknesses, except at the midplane, where the stresses result solely from the equilibrium 
contribution. There is no local mismatch in material properties at the midplane. As discussed 
earlier, the through-thickness extent of the mismatch contribution is controlled by the level of 
discretization used to model each layer and would decrease if the 0° and 90° layers were divided 
into sublayers. The most significant observation made from this figure is that the normal 
mismatch contribution is the component primarily responsible for the differences in the stress 
predictions provided by the present solution for the two laminates. The equilibrium and shear 
mismatch contributions have similar influences on the total distribution in both laminates, but 


the normal mismatch contribution has opposite effects. 




<*33 

(a) [0 4 /90 4 ] s Laminate 



<*33 


(b) [90 4 /0 4 ] s Laminate 

Figure 5.5. Contributions to o 33 at X 2 /b = 0.999 in [0 4 /90 4 ] s and [ 904 / 04 ], Laminates 
Extension 
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The influence of the normal mismatch effect on the total interlaminar normal stress is more 
vividly illustrated in Figure 5.6 which presents individual contributions to the a 33 stress 
component along the 0/90 interface of the [90 4 /0 4 ] s laminate. The mismatch effects are not as 
significant for the [0 4 /90 4 ] s laminate and are not discussed here. Distributions for this laminate 
are provided in Appendix C (Figure C.l) for completeness. Recall, the shear mismatch 
contribution to the normal stress is zero at all interfaces and therefore is not shown in the figure. 
Near the free edge, the normal mismatch component has magnitude approximately equal to the 
equilibrium contribution but of opposite sign; the mismatch contribution totally changes the 
character of the interfacial stress distribution and causes a reversal in stress near the free edge, 
tending toward positive stress as the free edge is approached. This type of behavior cannot be 
predicted by the KL solution, because of the constraint imposed by the their solution that 
stresses in all layers decay at the same rate. At the midplane, the mismatch terms are zero, and 
the stresses result solely from equilibrium. The current theory, the finite element method and the 
KL solution all predict the normal stress to be distributed over the interfacial plane X 3 = 0 
(midplane) as shown in Figure 5.7. The KL solution requires the through width stresses at the 
0/90 interface to have the same form, as displayed in Figure 5.2b. 

Similar distnbutions for the interlaminar shear stress 023 are provided in Appendix C 
(Figures C.2-C.3). As was the case for the interlaminar notmal stress the noimal mismatch 
contribution is responsible for the asymmetry in the through-thickness distributions a 23 observed 
in Figure 5.3. The influence of the mismatch contribution on the shear stresses is directly 

evident from Figure 5.4, resulting in an intensification of the stress in both magnitude and 
gradient. 

Another issue bnefly addressed is that of solution accuracy as related to the number of 
eigenfunctions employed in the stress expressions. Several solutions have been developed 

a-% 









TABLE 5.4. Complementary Energy for Different Orders of Approximation 


Complementary Energy for Different 
Orders of Approximation 

[0 4 /90 4 ], 

[ 9 O 4 /O 4 ], 

N 

(jj^XlO 3 (lb-in.) 

N 

^XIO 3 (lb-in.) 

n=2 (KL) 

4.9696 

n=2 (KL) 

4.4233 

2 

4.8185 

2 

4.3325 

6 

4.5743 

6 

4.0304 

10 

4.429022 

10 

3.711539 


throughout the course of this investigation, with each modification increasing the number of 
terms incorporated in the assumed stress states. Recall that the latest modification includes 
fourteen parameters in the stress expressions for the most general laminate configurations. As 
previously indicated, this number reduces to ten for cross-ply laminates. Figure 5.8 illustrates 
the variation of the interlaminar normal stress a 33 with X 3 as computed by solutions employing 
two parameters (2 term solution, n=2), six parameters (n=6), and the full ten term solution 
(n=10). The KL solution and finite element solutions arc also presented for reference. Stress 
expressions for the two term solution are provided in Appendix B. The stress expressions for the 
six term solution are identical to those given in equations (3.89-3.90) without the last terms 
(terms multiplying A 3 ). As the figure shows, the predictions of the two term solution are nearly 
coincident with the KL solution. The six term solution provides improved results, relative to the 
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finite element predictions, but still hasn’t captured the differences in response of the 0/90 and 
90/0 laminates. The ten term solution, as previously mentioned, captures this asymmetry. 

Laminate complementary energy values have also been evaluated for the various degrees of 
approximation. These values are provided in Table 5.4. As the table shows, the n=10 solution 
gives the lowest energy of the solutions considered. No statement can be made at this point 
regarding solution convergence. 

5.1.2 Symmetric Laminates - Bending Load 

In this section distributions obtained from the present methodology and the KL solution for 
[0 4 /90 4 ] s and [904/0 4 ] s laminates subjected to uniform bending loads are presented. Through- 
thickness and interlaminar stress distributions are provided. As was done in the previous 
section, through-thickness plots are presented only for the top half of the laminate, but the stress 
symmetry conditions about the laminate midplane are different for bending than they were for 
extension. In the case of uniform bending, interlaminar normal stress is antisymmetric about the 
midplane, and interlaminar shear stress is symmetric. Interfacial distributions are again plotted 
as a function of the normalized distance X 2 /b from the laminate center, where b is the laminate 
half width. Unless otherwise noted, all stress components are normalized by MuH/21, where 
Mn ' s the far field moment, given from classical lamination theory, H is the laminate half 
thickness, and I is the moment of inertia of a unit width section of the laminate. The normalized 
stresses are denoted by o 33 and o 23 . 

Figures 5.9 and 5.10 show comparisons between the present solution and KL predictions for 
through-thickness distributions of the interlaminar normal and shear stresses at two locations 
near the laminate free edge. The most significant observation to be made from these figures is 
that both solutions predict more severe interlaminar normal and shear stresses in the [90 4 /0 4 ], 
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(a) [0 4 /90 4 ] s Laminate 


X 3 /h 



(b) [90 4 /0 4 ] s Laminate 


Figure 5.9. a 33 Near Free Edge in [0 4 /90 4 ], and [90 4 /0 4 ], Laminates in Bending 
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(a) [0 4 /90 4 ] s Laminate 


X 3 /h 



- 0.030 - 0.020 - 0.010 0.000 

<*23 


0.010 


0.020 


(b) [90 4 /0 4 ] s Laminate 


Figure 5.10. o 23 Near Free Edge in [0 4 /90 4 ] s and [90 4 /0 4 ] s Laminates in in Bending 
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laminate than in the [0 4 /90 4 ] s laminate for the same loading, with the maximum stresses in the 
90/0 configuration being nearly twice as large as those in the 0/90 layup. Physically this is 
expected and is fortunate, since the 0/90 laminate is much stiffer in bending and therefore 
preferred for such applications. The forms of the through-thickness distributions obtained by the 
two approaches are however different. The most noticeable differences are in the interlaminar 
normal stress predictions for the [90 4 A)4] S laminate at X 2 /b = 0.999, as shown in Figure 5.9b. 
The present solution predicts much larger through thickness gradients than the KL solution, but 
the KL results show interlaminar normal stress twice as large as the present theory at the 90/0 
interface. Much smaller differences are observed in the predictions of the two methods for the 
interlaminar normal stress in the 0/90 laminate at the same X 2 location (Figure 5.9a). In both 
laminates, the solutions are in much closer agreement slightly away from the edge at 
X 2 /b = 0.995. The interlaminar shear stress predictions of the two methods, on the other hand, 
still have different forms at X 2 /b = 0.990, as displayed in Figure 5.10. 

The distributions presented in Figures 5.9 and 5.10 suggest that the mismatch effects have a 
stronger influence on the stress distributions in the [90 4 /0 4 ], laminate than in the [0 4 /90 4 ] s 
laminate. This is more clearly illustrated in Figures 5.11 and 5.12 which show width 
distributions of both interlaminar stresses at the 0/90 interface. Figure 5.1 1 shows the interfacial 
distribution of 033 for both laminates. The two methods predict essentially the same results in 
the [0 4 /90 4 ] 5 laminate. However, the results are quite different for the [90 4 /0 4 ] f laminate. The 
KL solution predicts a much steeper stress gradient, a maximum stress approximately three 
times larger than the present theory, and a larger stress reversal away from the free edge. The 
mismatch component reduces the stress at the free edge and flattens the distribution as the 
distance from the free edge becomes larger. Differences are also exhibited in the shear stress 023 
(Figure 5.12). For the [90 4 /0 4 ], laminate the current solution predicts maximum stress twice as 
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large as KL. 

The above discussion in conjunction with results obtained for the extension case suggest that 
the additional mismatch terms included in the present theory, over those included in the KL 
formulation, are required for accurate analysis in some cross-ply laminates subjected to bending 
loads. Additional analyses and comparison with other solutions and experimental results, as 
they become available, are necessary to support or refute this claim. 

The contributions of the various physical effects for this problem are similar to the extension 
case and are not discussed. Total through-thickness and interfacial distributions of stresses along 
with the individual contributions of global and local mismatch components are provided in 
Appendix C (Figures C.4 - C.7) 

It is also of interest to compare the stresses developed in laminates when subjected to 
extension load with those that develop when the laminate is subjected to uniform bending. Here, 
a comparison is made between the stresses at the 0/90 interface in the [ 04 / 904 ], laminate. In 
order to do so, a slightly different normalization scheme is used for the bending stresses than was 
used in the previous figures. In the following figures the bending stresses at the 0/90 interface 
are normalized by the average of the longitudinal stress d M in the 0° and 90° plies obtained from 
classical lamination theory at that interface. A similar normalization was used for the axial 
extension load case and provides an indication of the severity of the interlaminar stresses relative 
to the in-plane far field stress 6 n . 

Figure 5.13 shows comparisons of the 0/90 interface interlaminar normal and shear stresses 
for the two load cases. The shear stress distribution is very similar for both load conditions, but 
the character of the normal stress distribution is somewhat changed. In the bending analysis, a 
larger maximum normal stress is noted, and a larger stress gradient shifted toward the free edge 



-0.005 


Extension 



(b) Interlaminar Shear Stress 

Figure 5.13. Comparison of Bending and Extension Stresses 033 and a 2 3 at 0/90 Interface for 
[O4/9O4 ] s Laminates 
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is observed. Further, the normal stress does not experience as large of a reversal away from the 
free edge in the extension case as in the bending case. The boundary layer width is 
approximately the same for the two load conditions. The maximum magnitude of both the 
interlaminar shear stress and the interlaminar normal stress occurs for the bending case. 

5.1.3 Unsymmetric Laminates - Extensional Load 

This section discusses results predicted by the present theory for extensional loading (N u ) of 
an unsymmetric [0 4 /90 4 ] t laminate. Through-thickness distributions of normalized classical 
lamination theory in-plane stress a 22 , for extensional and bending loads (to be discussed in the 
next section) are provided in Figure 5.14. The stresses for extensional loading are normalized by 
N n /2H, and the stresses for bending loads are normalized by M n H/2I. The linear variation of 
the in-plane stresses through the thickness for the extension loading results from a positive Kn 
curvature which develops in the plate because of the laminate-membrane flexural coupling. For 
the bending load, the laminate develops both curvatures ten and K 22 . 

Through-thickness distributions (for the entire laminate thickness) of the normalized stresses 
a 33 and a 23 near the free edge are shown in Figure 5.15. The bottom surface of the laminate is 
denoted by X 3 /h = -l and the top surface corresponds to X 3 /h = 1. As the figure shows, shear 
stresses are nearly symmetric about the 90/0 interface and the normal stress is close to 
antisymmetric. Changing the stacking sequence from [0 4 /90 4 ] t to [90 4 /0 4 ] t simply results in a 
change in sign of the shear stress distribution. The maximum tensile normal stress occurs at 
X 3 /H = 0.35, and the maximum shear stress develops at the 90/0 interface. 

Figure 5.16 illustrates the character of the intcrfacial distributions of normal and shear stress 
along the 0/90 interface. Note that a solution based entirely upon overall equilibrium, predicts 
zero interlaminar normal stress at this interface. Recall that the intcrfacial normal stress, or the 




Figure 5.14. Through-Thickness 622 for [0 4 /90 4 ] t Laminate for Bending and Extension Loads 
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(a) C33 at X 2 /b = 0.999 
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(b) CTaa at X 2 /b = 0.993 


Figure 5.15. a 33 and O 23 Near Free Edge for [O4/9O4], Laminate - Extension Load 
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normal stress at any z-Iocation predicted by the equilibrium solution, is proportional to the 
moment of the in-planc transverse stresses 6 22 about an axis parallel to the longitudinal axis X, , 
lying in the X]-X 3 plane of interest. This moment is represented by the term B 7 (k) in equation 
(3.90). As can be seen from Figure 5.14, moment equilibrium is satisfied by the lamination 
theory stresses at the 0/90 interface; that is B 7 (k) is zero there. Interlaminar normal stresses are 
therefore not required at this interface for sublaminate equilibrium, but arise solely from local 
mismatch effects. As it turns out, the mismatch contribution is not critical for the load case 
considered, but it may be significant for compression loading, or in a combined loading 
situation. For the interlaminar shear stress (Figure 5.16b), mismatch and equilibrium effects 
both contribute at the 0/90 interface. The mismatch contribution shifts the peak total shear stress 
toward the free edge and results in larger maximum stress and steeper gradient than predicted by 
a solution based on equilibrium considerations alone. As a result, shear stresses predicted by the 
present theory are distributed over a smaller portion of the 0/90 interface than an equilibrium 
solution would indicate. 

5.1.4 Unsymmetric Laminates - Bending Load 

In order to illustrate the influence of load conditions on the interlaminar stress state, an 
analysis has also been conducted for uniform bending of the unsymmetric [0 4 /90 4 ] t laminate. 
Comparison of through-thickness distributions for the extensional (Figure 5.15) and bending 
(Figure 5.17) load cases indicates larger through-thickness gradients and larger normalized 
interlaminar stresses for extensional loading. However, for the bending load case, the maximum 
normalized interlaminar normal stress occurs near the 0/90 interface, and is larger at the 0/90 
interface than in the extension case. This stress is tensile if the stacking sequence is reversed, or 
equivalently, if the direction of the applied moment is changed. Consequently, bending may 
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represent the more severe loading case, in terms of delamination, for this laminate configuration. 

Interfacial distributions for the contributions of the various physical effects are clearly 
delineated in Figure 5.18. In contrast to the extensional loading case, mismatch and equilibrium 
effects contribute to the normal stress at the 0/90 interface, but the equilibrium contribution to 
the shear stress is zero. The interlaminar shear stress at this interface results from mismatch 
considerations alone, and is distributed over the laminate width with the self-equilibrating form 
shown in Figure 5.18b. The mismatch contribution to the normal stress tends to flatten the total 
distribution and reduces the maximum stress at the free edge from that predicted by the 
equilibrium contribution. 

5.1.5 Unsymmetric Laminates - Combined Load 

The results of the previous two sections can be superposed to obtain results for a variety of 
combined load conditions, and examination of the figures together provides some insight into the 
types and combinations of loading that will magnify the interlaminar stresses or make them less 
severe. One loading condition of interest is that which produces constant strain en and zero 
curvature ten in the laminate, since this type of condition allows for direct comparison with 
symmetric laminates subjected to uniform end extension. Under such a loading, the laminate 
will assume a pure cylindrical shape, with curvature K 22 - As Norwood* 761 has discussed, the 
upper half of a [0/90], laminate is equivalent to an unsymmetric [0/90], laminate constrained 
from deflecting out of plane by the restraint u 3 (XiX 2 , 0 ) = 0 applied to its lower surface when 
uniformly extended. Comparison of the interlaminar stress distributions for these two laminates 
then gives an indication, within the limits of linear theory, of the influence of the out-of-planc 
deflections on the severity of the interlaminar stresses. 
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Stress distributions along the 0/90 interface of [0 4 /90 4 ] g and [0 4 /90 4 ], laminates are provided 
in Figure 5.19. The maximum interlaminar normal stress at the free edge is tensile in both 
laminates (Figure 5.19a). The distribution for the unsymmetric laminate is, however, flatter 
away from the free edge than in the symmetric laminate, and attains a lower peak value. The 
maximum o 33 predicted by the present theory in the unsymmetric laminate is 24% smaller than 
that obtained for the symmetric laminate. Similar observations are made for the interlaminar 
shear stress o 23 (Figure 5. 19b). In both laminates the shear stress rises sharply as the free edge is 
approached and then decreases rapidly to satisfy the traction free boundary condition. The stress 
gradient in the unsymmetric laminate is seen to be slightly larger than in the symmetric 
laminate, but as was the case for the normal stress, the maximum shear stress is smaller in the 
unsymmetric laminate. The largest a 23 predicted by the present theory is approximately 12% 
percent smaller in the unsymmetric laminate than in the symmetric laminate. The predicted 
trends are in agreement with those predicted by Norwood,* 76 ^ using a finite element analysis, and 
with his conclusion that interlaminar normal stress and shear stress "in the unsymmetric 
[0 4 /90 4 ] t laminate are relieved by out of plane deflections"^. 

5.2 Angle-Ply Laminates 

In contrast to cross-ply laminates, angle-ply laminates, with individual laminae oriented at 
angles +6 and -6 to the global axis isolate the influence of local mismatch in coefficient of 
mutual influence on the interlaminar stresses, since Poisson’s ratios arc identical for the +0 and 
-6 layers. Results are presented only for the extension of symmetric angle-ply laminates 
because of difficulties encountered when analyzing some angle-ply laminates subjected to 
bending load, and unsymmetric angle-ply laminates subjected to extension loads. Before 
presenting the results for the symmetric laminates a brief discussion of the nature of this problem 
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is provided. 

In laminates where there is large local mismatch in coefficient of mutual influence between 
adjacent layers, in combination with a curvature k 12 the numerical solution of the non-linear 
system of equations became unstable, and jumped back and forth between two widely spaced 
solutions. One of the solutions made sense physically, but the other led to extremely small 
decay rates. This was found to be the case for both symmetric laminates subjected to bending 
loads and unsymmetric laminates under uniform extension. Incidentally, this problem was not 
restricted to angle-ply laminates. The solution for a [0/30], laminate subjected to uniform 
bending, for instance would not converge, while no problems were encountered with a [0/75], 
laminate. The difference in these laminates is the magnitude of the local mismatch in coefficient 
of mutual influence and the curvature k 12 . 

Unfortunately, the source of the problem is still unknown, but a few possible explanations 
are proposed. First, it is possible that some inconsistencies may have been introduced into the 
analysis by the manner in which the end conditions are being applied. Currently, displacements 
in the axial direction (ui) are prescribed and the other displacements are taken to be zero, 
simulating the conditions that would be present in an end gripped specimen. This is in effect 
introducing an axial dependence on the stresses that might be more severe for cases where tc J2 is 
large. Another approximation made, that may loose validity when large curvatures are present, 
is the assumption that only half of the laminate width needs to be modeled. However, the KL 
solution runs for all of these laminates. This suggests that if the above mentioned 
approximations are the source of the problem, the errors introduced by making them are 
magnified by the additional terms associated with coefficient of mutual influence mismatch that 
are incorporated in the current solution. Or there may be an error in the computer program 
associated with terms involving the product of the coefficient of mutual influence mismatch and 
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linearly varying stresses due to K n . Finally, there is the possibility that the problem is purely 
numerical. 

5.2.1 Symmetric Laminates - Uniform Extension 

For the uniform extension of symmetric angle-ply laminates, the transverse in-plane stresses 
O 22 from classical lamination theory are zero. Thus, interlaminar normal stress O33 and shear 
stress 023 316 n °t required for equilibrium, and have been shown to be small compared with O 13 . 
The present theory and KL solution predict identically zero O23 and O33, so the discussion that 
follows will focus on the effects of fiber orientation and stacking sequence on interlaminar shear 
stress 013 . Results are presented for two different stacking sequences - clustered [+ 02 /— © 2 1* a °d 
alternating [(±0 ) 2 ] s - of angle-ply laminates with 0 = 10°. To achieve the same level of local 
mismatch contribution for all laminates, a through-thickness decay length of one ply thickness ( 
h = 0.005) was used. For the clustered laminate this decay length corresponds to representing 
each layer as two sublayers of equal thickness. The decay lengths and remaining constants in the 
assumed stress states are provided in Table 5.5. For the angle-ply laminates, the only constants 
of interest are 4>i» <t> 3 . ^ 2 , and Aj, and the results of the study indicate <|>i = <(> 3 . As previously 
discussed, this reduced number of constants is a reflection of the fact that the approximate 
solution predicts O 13 to be the only non-zero interlaminar stress component for angle-ply 
laminates subjected to extensional load. Classical lamination theory stresses for the two 
laminates are provided in Table C.3. 

Through-thickness distributions of the interlaminar shear stress O 13 at X 2 /b = 0.999 for the 
two different stacking sequences are provided in Figure 5.20. In this figure, and subsequent 
figures, stresses are normalized by the average applied far field stress, Nn/2H. Normalized 
stresses are denoted with an overbar. The results in Figure 5.20a for a [(+ 1 0 > 2 /(— 1 0)2 ]« show 
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TABLE 5.5. Solution Parameters for Angle-Ply Laminates 


Solution Parameters For Angle-Ply Laminates 

Constant 

[(±10)2], 

[IO2/— IO2I3 

*2 d/in) 

120.7 

109.3 

^2 

4.135 

3.727 

A,X10 -5 

2.145 

2.060 

h (in) 

0.020 

0.020 


reasonably good correlation between all three methods. Similar results were obtained for other 
"clustered" angle-ply laminates with stacking sequence of the type [(+ 0 ) 2 /(- 0 ) 2 ]s- Recall that 
the interfacial shear stress O 13 predicted by the KL solution is proportional to the through- 
thickness integral of the in-plane shear stresses 0 ^ above or below the interface of interest. This 
integral is represented by the force sum in equation (3.87). For the clustered family of 
laminates this force sum is non-zero at all locations X 3 , and sub-laminate equilibrium requires 
the interlaminar shear stresses be non-zero throughout the laminate thickness. 

On the other hand, in "alternating" angle-ply laminates (Figure 5.20b) equilibrium is satisfied 
at the second interface by the lamination theory stresses, that is is zero there. Consequently, 
interlaminar stresses arc not required for equilibrium and the KL solution predicts identically 
zero shear stress O 13 at this interface, while the modified solution of the present study and finite 
element model predict shear stress of considerable magnitude. These large local interlaminar 
shear stresses can be attributed to the large mismatch in r) xy „ at the -0/0 interface. Also note 
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(a) Clustered Layers 



(b) Alternating Layers 


Figure 5.20. On at X 2 /b = 0.999 for [(+10) 2 /(— 10>2 ], and [(±10) 2 1, Laminates - Extension 
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that Oj3 is large at the 10/- 10 interface, having a magnitude of 0.25 times the average applied 
stress. The stress is large because the the mismatch in t| xy x is large. As shown subsequently, the 
shear stresses decrease for larger fiber angles as the mismatch in T| xy , decreases. 

The relative contributions of the two physical effects for this problem are depicted in Figure 
5.21 where the total through-thickness distributions of <J|3 are presented along with distributions 
of the global and local mismatch contributions to the stress field. It is clear from Figure 5.21 
that for the layer discretization used to obtain these results, i.e. representing each layer as two 
sublayers with thicknesses equal to one ply thickness, there is no local mismatch contribution in 
the top and bottom halves of the first and second layers, respectively, whereas global mismatch 
or equilibrium contributes throughout. As discussed previously, the through-thickness extent of 
the local mismatch contribution to the stress field is governed by the thickness of the sublayers 
adjacent to interfaces where a mismatch in material properties is present. In contrast to the 
clustered configuration, both equilibrium and mismatch contribute throughout all layers of the 
[(±10)2], laminate (Figure 5.21b). Again, the extent of the mismatch contribution would be 
changed if different "sublayer" thicknesses were used. 

Interlaminar distributions of 013 along interfaces in the [(±10)2!* laminate are presented in 
Figures 5.22 and 5.23. Figure 5.22 shows the predictions of the present solution, the KL 
solution and finite element analysis along the first and second interfaces of the laminate. The 
three methods compare quite well at the first interface (Figure 5.22a), but only the present 
solution and finite element results agree along the second interface. As discussed above, the KL 
solution predicts zero shear stress O13 along the entire second interface, and the total stress 
results solely from local mismatch effects. In Figure 5.23 the distributions of the total shear 
stress and the individual mismatch contributions along the first interface are delineated. The 
local mismatch distribution is self-equilibrating, integrating to zero over X2, while the 




(b) Alternating Layers 

Figure 5.21. Contributions to c i3 in [(+ 1 0) 2 /(— 1 0) 2 ], and [(±10) 2 ], Laminates - Extension 
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equilibrium contribution is equivalent to a non-zero intcrfacial force. The results indicate a 
larger contribution from local mismatch at the intersection of the interface and free edge than 
from global mismatch (equilibrium). However, the boundary layer width is approximately the 
same for both effects. 

The influence of fiber orientation on the stress field characteristics is illustrated in Figures 
5.24-5.26 for both "clustered" and "alternating" stacking sequences. Figure 5.24 shows the 
variation in coefficient of mutual influence and Poisson’s ratio as a function of 0 for a T300- 
5208 graphite epoxy with the material properties given in Table 5.1. As can be seen, T| 12 . x 
attains a maximum value at 0 = 15°. Thus, the interlaminar shear stresses are expected to be 
maximum in angle ply laminates with adjacent +15/-15 layers. Figure 5.25, which shows the 
variation in the maximum intensity of O 13 with 0 , for both stacking sequences, demonstrates that 
this is in fact the case. The parameters in the assumed stress expressions are also functions of 0. 
The variation of 4>i with 0 for the two stacking configurations is depicted in Figure 5.26a where 
it is evident that <f»i attains a minimum at 0 = 30° for both stacking sequences. The constant A] 
on the other hand reaches a maximum at 0 = 25° (Figure 5.26b). The constant % 2 . is independent 
of 0,but does depend on stacking sequence. As indicated in Table 5.5 =4.135 for the 

alternating sequence and X 2 = 3.727 for the clustered sequence. 

5.3 Quasi-Isotropic Laminates 

Quasi-isotropic laminates are currently used extensively in practice because of their isotropic 
in-plane elastic properties. Symmetric quasi-isotropic laminates made with equal percentages of 
0, 90, +45, and -45 degree laminae orientations are examined here. There are twelve unique 
configurations of this type, as shown in Table 5.6, if it is assumed that the +45 and -45 layers can 
be interchanged. In the table, the laminates are divided into two groups; those with adjacent ±45 
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layers and those with interspersed ±45 layers. This terminology will be used throughout the 
subsequent discussion. 

In quasi-isotropic laminates, in contrast to the special cases of angle-ply laminates and 
cross-ply laminates discussed in the previous sections, both a mismatch in v 12 and r|i 2 ,i may 
exist between adjacent layers. Consequently, all stresses are in general non-zero, and the total 
fourteen parameters <J>i(i = l,5), (J = 1,6), and A k (k= 1,2,3), in the assumed stress 
expressions must be determined in order to evaluate the stresses in the individual plies. 

5.3.1 Extensional Load 

Analyses have been conducted for the uniform extension of all twelve laminates listed in 
Table 5.6, but because of the large number of plots necessary to characterize the stress field near 
the free edge in each laminate, results in the form of through-thickness and interfacial stress 
distributions are provided for only three of these: [90/45/0/—45],, [45/90/0/— 45], and 
[0/45/-45/90],. These laminates were chosen to illustrate the influence of adjacent ±45 layers 
and interspersed ±45 layers on the interlaminar stress state and to study the relative magnitudes 
of the mismatch and equilibrium contributions to the stress field for different stacking sequences. 

Stress distributions for the interspersed [90/45/0/-45], and [45/90/0/-45], laminates and the 
[0/45/— 45/90], laminate with adjacent ±45 layers are provided in Figures 5.27-5.31. Stresses 
determined by finite element analysis, the KL solution, and the present theory (equations (3.86- 
3.89)) are shown. The finite element results for the [90/45/0/— 45], and [0/45/-45/90], laminates 
were generated previously by Herakovich [9] using the program CLFE2D, and the finite element 
results for the [45/90/0/ — 45] a were obtained by the author using Norwood’s 1761 program. The 
present theory results were generated using a through-thickness decay length of one layer 
thickness. All layers were of equal thickness, t w = h = 0.005 in. The parameters ft , Xj, and A k 
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TABLE 5.6. Quasi-Isotropic Laminates 


Interspersed ±45 Laminates 

Adjacent ±45 Laminates (ksi) 

[90/45/0/— 45], 

[90/0/±45], 

[0/ — 45/90/45 ], 

[90/±45/0], 

[45/90/0/— 45], 

[0/90/±45], 

[45/90/— 45/0], 

[0/±45/90], 

[45/0/90/— 45], 

[±45/90/0], 

[45/0/— 45/90], 

[±45/0/90], 


for the laminates considered in this section are provided in Table 5.7. Solution parameters for 
the remaining quasi-isotropic laminates listed in Table 5.6 are provided in Appendix D (Tables 
D.2-D.5). In all figures, stresses are normalized by the average far field stress On =Nn/2H 
where Nji is the far field load, and 2H is the laminate thickness. Classical lamination theory 
stresses for the three laminates subjected to uniform extension are provided in Table D.l. 

Figures 5.27-5.29 display through-thickness distributions of the interlaminar stress 
components for the two interspersed stacking sequences. The results in Figure 5.27 and 5.28 for 
the interlaminar stresses in the [45/90/0/-45], laminate show similar trends in the predictions of 
all three methods. The results of the present theory for the interlaminar shear stress 013 , shown 
in Figure 5.27a, and the interlaminar shear stress O 23 , presented in Figure 5.28, however, show 
better correlation with the finite element solution than exhibited by the KL solution. The present 
solution more accurately predicts stress magnitudes, and more closely approximates the 
through-thickness variations in the stresses than docs the KL solution. 
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Similar observations are made with regard to the normal stress predictions, depicted in 
Figure 5.27b. All methods correlate fairly well, except at the laminate midplane. At the 
midplane, the present solution and KL solution predict larger compressive interlaminar normal 
stress than finite elements. A possible cause for the difference between the finite element 
predictions and the present solution is discussed later in this section. 

The similarities in the predictions of the three methods for stresses in the [45/90/0/-45], 
laminate discussed above suggests that, for this laminate, mismatch effects do not have much of 
an effect on the overall distribution of stress in the laminate; that is, global equilibrium 
dominates the development of the interlaminar stresses. 

In the remaining interspersed laminates, where the 90 degree laminae is positioned at the top 
or bottom of the stack, [90/45/0/— 1 45], and [45/0/— 45/90],, mismatch effects have a significant 
influence on the total stress distributions. This is clearly illustrated for the [90/45/0/— 45], 
laminate in Figures 5.29 and 5.30. These figures show comparisons of the three methods for 
through-thickness distributions of 0| 3 (Figure 5.29a) and interlaminar normal stress (Figure 
5.29b), and for interfacial distributions of O 23 (Figure 5.30). Evident from Figures 5.29a and 
5.29b is the significant improvement of the present theoiy over the KL solution, with the present 
theory predicting trends more similar to those displayed by the finite element results. The 
modified solution of the present study and the finite element model predict more severe 
through-thickness gradients for the interlaminar normal stress, and generally larger shear stresses 
013 than the KL solution. At the first interface (90/45), for instance, the present theory and finite 
elements predict interlaminar shear stress O13 (Figure 5.29b) to be approximately 6% of the 
average far field stress, where the KL solution predicts identically zero stress. At this interface, 
equilibrium is satisfied by the KL solution, due to the in-plane shear stress O 12 being zero in the 
90 degree layer. The non-zero stress 013 at this interface then results solely from the mismatch 
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Figure 5.27. 0,3 and o 33 at X 2 /b = 0.999 in [45/90/0M5], Laminate - Extension 
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in "H 12. 1 between the 90 degree and 45 degree layers. 

There are also differences in the predictions of the three methods for the interlaminar shear 
stress 023- Through-thickness finite element data is not available for but interfacial results 
(first and second interfaces), presented in Figure 5.30, indicate that the present solution compares 
more favorably with finite element results than do the KL results. Recall, that the KL solution 
requires the stress distributions at these interfaces to be of the same form, since stresses arc 
constrained to decay at the same rate in all plies. The present theory and the finite element 
method, on the other hand, exhibit variations in the character of the interfacial distributions for 
these two interfaces. The most significant difference in the interfacial predictions for O23 
observed is at the second interface where the local mismatch contribution causes a reversal in 
stress near the free edge (Figure 5.30b). These differences and the differences discussed above 
in the o 13 and a 33 predictions could be significant when evaluating the structural response of 
candidate laminates to an applied load, and are clear evidence of the influence of local mismatch 
effects on the total stresses and the necessity for including these effects in the stress analysis. 

Similar conclusions can be drawn for the quasi-isotropic laminates with adjacent ±45 layers; 
that is, in these laminates local mismatch effects are significant and a design or analysis based 
upon stress predictions obtained from global equilibrium considerations alone can be 
misleading. This observation is illustrated in Figure 5.31, which shows predictions for 
interlaminar shear stress Oj 3 and a 33 in a [0/45/ — 45/90] a laminate which arc typical of results 
obtained for laminates having adjacent ±45 degree layers. Note from Figure 5.31a that the 
maximum shear stress develops at the 45/-45 interface, where the mismatch in r| 12>1 is the 
largest. The present theory predicts a 13 at the 45/-4S interface 5% smaller than predicted by 
finite elements, while KL underpredicts the finite element stresses by 38%. Similar trends were 
observed in the predictions of the three methods for the other quasi-isotropic laminates with 
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adjacent ±45 layers. In the other laminates, the shear stress predicted at the +4S/-45 interface by 
the present solution and the finite element method was as much as 50% to 100 % larger than 
predicted by the KL solution. 

Results are also provided in Figure 5.31b for the interlaminar normal stress in the 
[0/+45/-45/90], laminate. Again, these results are typical of those obtained for the group of 
adjacent quasi-isotropic laminates. Note that the finite element solution predicts larger through- 
thickness gradients in the normal stress than predicted by cither the KL solution or the present 
solution. The present theory appears to be starting to pick up the trends displayed by the finite 
clement results, but does not predict the large compressive stress at the 45/-45 interface. 

One possible cause for the differences in normal stress predictions of the present theory and 
finite element methods at the 45/-45 interface of the laminate shown in Figure 5.31b is the 
assumed lack of coupling, in the present analysis, between the coefficient of mutual influence 
mismatch and the interlaminar stresses a 23 and o 33 , and between Poisson’s ratio mismatch and 
interlaminar shear stress o ]3 . This lack of coupling is also thought to be the reason for the poor 
correlation in the normal stress predictions of the present solution and the finite element results 
at the midplane of the [45/90/0/— 45] s laminate discussed at the beginning of this section, and for 
the differences in the normal stress predictions of the present theory and finite elements at the 
first interface of the [45/-45/0/90], discussed in Section 4.1.3 (Figure 4.6). 

Recall from the development in Chapter 3, that the coefficient of mutual influence mismatch 
was assumed to affect only the Oj 2 and Oj 3 components of stress while the Poisson’s ratio 
mismatch was assumed to affect only the o 22 , 023 and a 33 stress components. This assumption 
could be made for the problem studied since the assumption that stresses are independent of the 
longitudinal coordinate Xj led to a reduced system of equilibrium equations, where the stresses 
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Oi 2 and Oo uncoupled from O 22 . <*23 and <* 33 - However, as briefly discussed in Chapter 3, the 
stresses are coupled by the compatibility equations, and would be in the equilibrium equations 
for more general problems where the longitudinal independence could not be assumed. 
Therefore, an improved stress field assumption would include the same eigenfunctions in the 
expressions for all of the stresses, and would relate interlaminar stresses O 22 , <* 23 . and <*33 to the 
mismatch in coefficient of mutual influence and interlaminar shear stress O 13 to Poisson’s ratio 
mismatch. 

To summarize the results of this section, the above discussion indicates that stacking 
sequence has a significant influence on the relative magnitudes of the mismatch and equilibrium 
contributions to the stress field. In the laminates with interspersed ±45 layers, the factor having 
the largest effect on the magnitude of the different contributions is the location of the 90 degree 
laminae in the stacking sequence. This is explained by considering Poisson’s ratio mismatch. 
The variation in Poisson’s ratio with fiber orientation 0 for a T300-5208 graphite-epoxy, is 
shown in Figure 5.24. As can be seen from the figure, the mismatch in Poisson’s ratio between a 
445 degree or -45 degree layer and a 0 degree layer is much larger than the mismatch between a 
45 degree layer and a 90 degree layer. In the interspersed laminates with the 90 degree laminae 
as a middle layer, there is only one occurrence of adjacent 0 degree and 45 degree layers. In the 
interspersed laminates with the 90 degree laminae placed at the top or bottom of the stack, on the 
other hand, there are two occurrences of adjacent 45 degree and 0 degree layers, and the 
mismatch contribution to the total stress state is more significant. In the laminates with ±45 
layers, the coefficient of mutual influence mismatch is large, and the local mismatch effects 
again contribute significantly to the total stress. These observations are made clear by 
examining Figures 5.28, 5.29 and 5.31 together and comparing the relatives differences in the 
predictions of the present method and KL solution for the different stacking sequences. 
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53.2 Bending Load 

In order to illustrate the influence of load condition on the interlaminar stress, quasi-isotropic 
laminates subjected to bending load (negative curvature k m ) have also been analyzed. In this 
section, comparisons of the through-thickness distributions of the interlaminar normal stress and 
shear stresses for a [45/-45/90/0], laminate subjected to bending (-K n ) and extension loads are 
presented (Figures 5.32-5.34). As usual, distributions are presented for the top half of the 
laminate. In the case of bending, interlaminar normal stress is antisymmetric about the 
midplane, and interlaminar shear stress is symmetric. For the extensional load case, interlaminar 
shear stress is antisymmetric about the midplane and interlaminar normal stress is symmetric. 
As has been done in the previous sections, bending stresses are normalized by Mu H/2I, and 
extensional stresses are normalized with respect to Nu/2H where Mu and N n are the applied 
far field loads. All results were generated using a through thickness decay length of one layer 
thickness (h = 0.005 in.), and the solution parameters ft, Xj, and A k given in Table 5.8. 

Comparison of the through-thickness stress plots in Figures 5.32-5.34 shows similar 
distributions for all interlaminar stress components for bending and extension load. In both 
cases, the maximum interlaminar normal stress is tensile and develops at X 3 /h = 1.6. The 
maximum shear stresses also occur at the same through-thickness location in both load cases. 
The maximum interlaminar shear stress o 13 develops at the first interface (X 3 /h = 3), and the 
interlaminar shear stress attains a maximum at the second interface. The combined load 
case would then represent a more severe condition in terms of dclamination potential, with 
bending and extension stresses combining above the laminate midplane. 



Solution Parameters For [45/-45/90/0], Laminate 


Constant 

Extension 

Bending 

<h (1/in) 

57.07 

72.03 

<t >2 (1/in) 

63.25 

58.37 

<t>3 0/in) 

99.07 

69.22 

<i>4 d/in) 

66.82 

46.96 

<t>5 (1/in) 

95.58 

48.99 

X, 

1.966 

4.744 

^2 

4.034 

5.372 

X 3 

5.070 

10.71 

^4 

18.02 

11.80 

h 

1.565 

4.816 

*6 

6.062 

7.777 

A 2 X10“ 5 

(psi) 

1.773 

1.388 

A 2 X10~ 5 

(psi) 

-3.741 

-2.372 

A 3 X10 -5 

(psi) 

-1.892 

-2.857 
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Figure 5.32. o 13 at X 2 /b = 0.999 in [45/-45/90/0]* Laminate, Bending and Extension 



135 


X 3 /h 



Figure 5.33. 033 at X 2 /b = 0.999 in [45/— 45/90/0], Laminate, Bending and Extension 
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CHAPTER 6 

CONCLUSIONS AND RECOMMENDATIONS 
6.1 Summary of Method 

The purpose of this study was to develop an efficient approximate solution for interlaminar 
stresses near free edges of finite width, laminated composites subjected to bending and extension 
loads and to demonstrate the utility of the methodology. The analysis developed is an extension 
of a method recently presented by Kassapaglou and Lagace^ 20 ^ and is based upon the principle 
of minimum complementary energy and an assumed stress state, derived by considering material 
mismatch considerations and global equilibrium requirements. In the KL solution only the 
mismatch between laminae and laminate material properties was considered. The present 
solution extended their technique by including additional terms in the stress assumptions, which 
account for the local material property mismatch in coefficient of mutual influence and Poisson’s 
ratio between adjacent layers in the laminate. The differential equations of equilibrium, the 
interfacial traction continuity and boundary conditions of stress were identically satisfied by the 
assumed stresses. The strain compatibility equations and interfacial displacement continuity 
conditions were satisfied approximately by imposing the stationary condition of laminate 
complementary energy. 

The developed methodology is general, and in theory, can be extended for the analysis of 
other load cases, such as in-plane bending and torsion, or to more general structural 
configurations, provided a plane stress state is recovered in the member’s interior region, and the 
stress components do not vary with the longitudinal coordinate. The plane stress solution can be 
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obtained from an analytical solution, if available, or from a numerical technique, such as the 
finite element method. The methodology could then be employed in a global-local analysis to 
obtain stress solutions in regions of high stress gradients with a coarser solution used to define 
the response outside of these regions. 

As the previous paragraph indicates, a major advantage of the method presented is the fact 
that in-plane stresses, obtained from classical lamination theory, are the only required input to 
the solution. Other advantages, as compared with numerical solutions, or some of the 
complicated analytical models that have appeared in the literature, include the relative simplicity 
of the theory, in terms of the number of parameters that must be determined in order to obtain 
stress distributions, solution efficiency, output readability, and the ease of application for the 
simple geometry considered. For the most general laminate, only fourteen parameters must be 
determined prior to calculating stresses. This number reduces to ten for cross-ply laminates and 
to three for symmetric angle-ply laminates subjected to extensional load. For bending of angle- 
ply laminates, or extension of unsymmetric angle-ply laminates, the number reduces to six. In 
any case, the number of parameters is independent of the number of layers in the laminate, their 
material properties, and orientations. Consequently, the analysis can be applied to laminates 
with a large number of plies. Finally, a major asset of the solution is the insight it provides into 
the fundamental physical mechanisms, global equilibrium and local mismatch effects, that 
contribute to interlaminar stress development. Contributions from each of these effects are 
clearly delineated in the solution output. This type of information cannot be obtained from a 
finite element analysis. 

6.2 Conclusions 

The development of the solution methodology has been accompanied by application of the 
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stress analysis to several finite width laminates subjected to uniform extension and bending load. 
Symmetric and unsymmetric cross-ply laminates, symmetric angle-ply laminates, and symmetric 
quasi-isotropic laminates have been studied. Comparisons in stress predictions were presented 
for results obtained using the present solution, the KL solution, and finite element solutions, 
where available. Several conclusions can be drawn based upon the results of these studies: 

(1) The present method compares well with finite element methods and provides 
significantly improved stress predictions as compared with the KL solution which is based 
entirely upon global equilibrium considerations. In particular, the present solution more 
accurately predicts stress magnitudes and interlaminar stress gradients near interfaces where 
there is a large material property mismatch between adjacent layers. Further, the present method 
accurately characterizes the through-thickness gradients in the layer stress fields near the 
interfacial surfaces, while the KL solution predicts the laminate behavior to be qualitatively the 
same throughout the thickness of the layers. 

(2) The relative importance of local mismatch and equilibrium considerations was found to 
be a function of stacking sequence. This was the case for all of the laminate families analyzed. 
Specifically, the following observations were made: 

(a) In angle-ply laminates subjected to axial loading large mismatches in the 

coefficient of mutual influence induced large interlaminar shear stress In 

alternating stacking sequences, [(+©/— 0) 2 ] s . local mismatch resulted in stresses of 
considerable magnitude at interfaces where the KL solution predicted zero stresses. 
These results were in agreement with finite clement results. 

(b) In symmetric cross-ply laminates, subjected to both bending and extension 
loads, local mismatch effects had a more pronounced effect on the stress predictions in 
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[90/0] s laminates than in [0/90], laminates, and were the cause of the asymmetry in 
stresses observed in these laminates. More specifically, the mismatch contribution 
associated with Poisson’s ratio mismatch and the direct assumption on interlaminar 
normal stress (referred to as normal mismatch contribution) was found to be the 
component primarily responsible for differences in the stress predictions provided by the 
present solution for the two laminates. 

(c) In the quasi-isotropic laminates considered (laminates with equal percentages of 
0, 45, -45, and 90 degree layers) the degree of local mismatch contribution was also 
found to be dependent on stacking sequence. Specifically, in interspersed 
configurations, with the 90 degree layer separating the 45 and -45 degree layers, the 
mismatch contribution to all interlaminar stress components was small compared with 
the equilibrium contribution. For these laminates it was concluded that the magnitude of 
the stresses was primarily dependent on the force and moment developed at any 
through-thickness location by the intralaminar stresses O 22 . This was physically 
explained by the fact that the Poisson’s ratio mismatch and mismatch in coefficient of 
mutual influence between adjacent layers in this laminate are smaller than in 
interspersed laminates with the 45 and -45 degree layers separated by zero degree layers 
and in laminates with adjacent +45 and -45 degree layers. These results clearly 
indicated the influence of local mismatch effects on the stresses and the necessity for 
including mismatch considerations in laminate design. 

(3) For the laminates considered, bending and extension loads resulted in similar 
interlaminar stress distributions in terms of magnitude and boundary layers widths. 


(4) Comparison of results for an unsymmetric cross-ply [0/90], and a symmetric cross ply 
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[0/90], subjected to extension load showed that the out of plane deflections of the unsymmetric 
laminate relieved both interlaminar shear stress and interlaminar normal stress. 

(5) Difficulties were encountered in the analysis of laminates having large local mismatch in 
coefficient of mutual influence in combination with curvature k 12 . This was true for both 
symmetric laminates subjected to bending loads and unsymmetric laminates under uniform 
extension. Thus, at the present time, there is a limitation on the general applicability of the 
developed solution. The source of the problem has not been identified but may be attributed to 
some possible inconsistencies introduced by the manner in which end conditions are currently 
applied, or possibly with the assumption that only half of the laminate width needs to be 
analyzed. The assumption of modeling half of the laminate width is recognized to be an 
approximation for laminates with off-axis plies. This approximation may loose validity when 
large curvatures k 12 are developed. It is also possible that the problem is purely numerical or is 
the result of a coding error. This topic was discussed in more detail in Section 5.2 and is an area 
requiring additional study. 

6.3 Recommendations for Future Work 

The solution method developed has been shown to generally produce accurate predictions for 
interlaminar stresses near straight free edges of laminated plates. Although the model was 
developed for specific load conditions, and was based upon some assumptions that limit its 
application, it can be extended to more general analyses. Some recommendations for extension 
of the method and additional possible applications are presented in this section. First, however, 
some suggestions are provided for solution modifications that might lead to increased accuracy 
in the stress predictions. Recommendations related to solution efficiency are also provided. 


The results of the previous chapter showed that the least accurate results were obtained for 



142 


the interlaminar normal stress. It was suggested that improvements in the stress predictions 
might be achieved by incorporating additional exponential terms in the assumed stress 
expressions which account for the coupling between coefficient of mutual influence mismatch 
and interlaminar stresses CT23 and CT33, and the coupling between Poisson’s ratio and interlaminar 
shear stress O13. This could be accomplished by including an additional term in Equation (3.87) 
for the interlaminar shear stress O13 with the same form as that multiplying Aj but proportional 
to the mismatch in Poisson’s ratio. Similarly, additional terms in the expressions for O23 and 
033 with the same forms as the mismatch terms incorporated presently, but proportional to 
mismatch in coefficient of mutual influence could be included. The remaining stress 
components would then by determined from the differential equations of equilibrium. 

Another possible improvement to the stress assumptions would allow for a priori unspecified 
through-thickness decay rates for the local mismatch effects. One way in which this could be 
incorporated into the present theory would be to assume exponential functions for the through- 
thickness variations of the local mismatch, expressed in terms of unknown through-thickness 
decay parameters as was done for the width variations in the present solution. The O13 variation 
would require one exponential function, while 023 and O33 would require a combination of two 
functions. 

The above recommended modifications should lead to improved accuracy, however, the 
additional level of complexity introduced, may not be worth the effort. Further, solution 
efficiency will decline as the number of terms is increased. In this connection, the issue of 
computational efficiency of the current solution is addressed. The current solution generally 
executes in less than 15 CPU seconds on an IBM RS6000, for twenty different initial 
approximations to the solution. This run time, is of course dependent on the convergence 
tolerance used when solving the non-linear system and how good the initial guess is. In 
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generating the results presented in this thesis, as many as 500 initial guesses were employed for 
the most general laminates. This large number of starting points was felt necessary because of 
the difficulty with the present method of determining a good initial guess. Additional work is 
necessary to improve this aspect of the solution. 

Other areas of possible future work include extension to additional load cases such as in- 
planc bending and torsion, and to problems where stresses arc a function of the longitudinal 
coordinate X|. The method could also be extended to analyze more general structural 
configurations. One configuration of particular interest, is a stiffened panel. Extension of the 
method to analyze skin-stiffener interface stresses, or stresses in the tip of a stiffener blade is an 
area recommended for future work. 

Coupling of the stress analysis with an experimental program and some type of failure 
analysis or delamination initiation predictions is also recommended. Analytical failure 
predictions or delamination initiation predictions, based upon the stress analysis of the present 
work, could then be correlated with experimentally observed behavior. 

Finally, the methodology should be linked with a numerical optimization program to develop 
a design capability for laminated composites in which interlaminar stresses are considered. This 
capability will provide engineers with an efficient methodology for designing delamination 


resistant structures. 
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In this appendix the energy expression presented in Chapter 3 is expanded. To facilitate the 
expansion contracted notation rather than tensor notation is used for the stresses. The assumed 
stresses are written in the form 


f’mjGml + Fm2G m 2 + Pm3^m3 F m 4 G m 4 + F m 5 G m 5 


where F^ are functions of the decay parameters and y, and Gnu are functions of z. The 
functions G m i can be taken directly from the stress expressions in equations (3.86)-(3.90) and are 
not repeated. 


F ml =0 m=2,3,4 

T7__ _ 


F 5 i=a 5 ,e^ ,y 
^61 = (1 + a^iel) 


F 22 = ( 1 + a 2 2e" <l>2y + ^e ’^ y ) 
Fm2 = ( a m2 e"^ 2y + a^e - ^^) 

F m2 =0 

m = 3,4 
m = 5,6 

Fm3 =0 

m =2,3,4 
m = 5,6 

F„4 - 

F m4 =0 

m = 2,3,4 
m = 5,6 

F m 5 + a n ,ioe _Xs * sy + 

FmS =0 

m = 2,3,4 
m = 5,6 


where the a„ ’s are deteimined from inspection of equations (3.86)-(3.90). 


The energy expression can be written as given below 


n c = C[ + A] C2 + A2C3 + A3C4 + A] A2C5 

A2A3C6 + Ai A3C7 + AjC8 + A2C9 + A3C10 


where the terms Q are defined below 
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C, = \mh%72 + f 23 823 + f 24 8 24 + 4 861 + fllSll> 

+ ^22^32822832 + ^61^32861 832 + f42 f5 1 842851 + f6l(2286l822 
+ (a22/<t>2 + a 23 /( <t>2^2 ) )822 +( 361/4*1 )861 

Q = fsi ^53^51 853 + (si (63861 863 + ^63^22863822 + f63f32863832 + ^53^42 853 g42 

+ [ a64/4*3 + 365 /( ^-24*3 ) 1863 

C 3 = ^22^24822824 + f32f34832834 + f42(44 g42g44 + f22f34g22g34 

+ f24f32g24g32 + ^61 f24§6lg24 + (61 (34 861834 + (5lf44g5l844 

+ [a26/4>4 + a27/(^34>4) + 328/(^44*4) 1 g24 

C4 = f22f25g22g25 + (32(35 832 835 + (*2(45 842845 + (22(35822835 
+ (25 (32 625 832 + (51 (25 861 825 + (si(3586l835 + (5lf45g5lg45 
+ [ a29/4>5 + 3210/(^-54*5) + a 2l l/(^6 4>5) 1 825 

C5 = ^63^24 863 824 + (63(34863834 + f53f44g53844 

C6 = f24f35g24835 + f34f25 834 825 + f24f25824825 + f34f35834g35 + (44(45 g44g45 
C7 = (53(25863825 + (53f35 863 835 + f53f45853845 
C 8 = l/2[ 4g63 + f53g53 ] 

C 9 = l/2[ f 2 4 g 24 + (34 g34 + f44g44 ] + f24(34 g24 g34 
C 10 = l/2[fl 5 g2 5 + (35835 + (45845 ] + (25(35825835 


The integrals required to evaluate the energy expression arc expanded below. In these 
expressions lower case f m j and g m ; arc used to denote the integrated product o( the F m j and G m , 
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functions. 


form = n = 2,3,5: 


JjRmnO m a n dydz 


fm 2 Sm 2 2A2f m 2fm4gm2goi4 + 2A3f m 2f m 5goi2goi5 +2 A 2 A 3 f m 4 f m 5 
+ A2fm48m4 + Affj^gmsj Roin/2 


for m=2, n=3: 


JjRmnCJmOndydz 


^ ^22^32822832 + ^2 (f22 ^34 822 834 + f24f32g24g32) 

+ A 3 (f 2 2 f 35 g 22 g 35 + hshlZlSl&l) + A 2 A 3 ( f 24 f 35 g 24 g 35 + f 25 f 34 g 25 g 34 ) 
+ A 2 ( f24 ^34 824 834 ) + A 3 (f25 f 35 g25 835 )] ^mn 


form = n = 5,6: 


JJ^mn^m^n^ydz — fmlgml 2 Aj (f m i f m 3g m ig m3 ) + A| (f m3 g m 3^ Rj^/2 


for (m = 6,n = 2,3),(m = 5,n = 4): 


JJ^mn^m^n^ydz — ^ f m i f^goil gnj (^3^280138112) ^2 (^ml fn 4 8m 1 Sn 4 ) 

-h A 3 (f m i f^gmi g^) + Ai A 2 (fm 3 fn 4 goi 3 gn 4 ) + Aj A 3 (f m 3 fn 5 gm 3 gn 5 )j 


where 


f 22 = 2a2i/<t>2 + 2322 /(^ 1 ^ 2 ) + a 22 e 2,2 + 2a22 a 23 e 2,3 + a 23 e 3,3 
for (m = n = 3,4), (m = 2,n = 3) 

fm2fri2 = ^ 2 (^ 2 ^2,2 4* 3n3^3,2) a m3 (^2^2,3 ^3^3,3) 

^22^24 = a 26/$4 + a 27/(^3^4) + a 28/(^44>4) + a 22( a 26 e 2,6 + a 27 e 2,6 + a 28 e 2,8) 
+ a 23( a 26^3,6 + a 27^3,7 + 328^3,8) 
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for (m = n = 3,4), (m = 2,n = 3) 

fm 2 f n 4 = ^m2 (^^2,6 + ^nl^2,6 + 811862,8) + 81113(^62,6 + 8 n 763.7 + Sj^.g) 


h2hs - a 2 9/<l>5 + 8210/(^54*5) + 821 1/(^6 4*5) + 822(82962,9 + 821062,10 + 821162,11) 
+ 823 ( 82963,9 + 821063,10 + 821063 , 10 ) 

for (m = n = 3,4), (m = 2,n = 3) 

f m 2fn5 = a m 2(8n9C2,9 + 8nioC2,io + ani1e2.11) + a m 3(a n 9e2,9 4-a n ioe3,io + 8ni1e3.11) 


for (m = n = 2,3,4), (m = 2,n = 3) 

fm4fn5 = 8tn6(8t,9 e 6,9 + 8nl0 e 6,10 + 8 nlie 6 .ll) 4- am7(8n9 e 7,9 + 8„1066.10 + 8 nl 1 e 7 . 11 ) 

+ 8 m 8(a n 96g,9 + Snioegio + anll e 8.1l) 

f m 5fn4 = 8 n# (a m 9e 6t 9 + amioe6,io + amii e6,n ) + a n 7(a m 9e7 >9 + amioe6.io + ami 1 e 7 . 11 ) 

+ a^an^g^ + amioegio + ami^n) 


for (m = n= 2,3,4) 

(m4 fn4 8m6(ani666,6 4" 8m766,7 "1” 8m866,8) 4“ 8m7 (8m666,7 4” ^m7^7,7 4" 8 m ge7,g) 

4" a m g(a m6 e 6 , 8 + a m 7 e 7 ,8 + amgeg.g) 


for (m = n = 2,3,4) 

fm 5 fn 5 = a m 9(a m 9C9 9 + amioCg ,10 4 - 8m9e9.11) + a m io(a m 9C9,io 4 - amioCiojo 4 - amiiCio.il) 
4 - am 11 (8m9 69, 11 + aml0 e 10,ll 4 ^ 8mll e ll,ll) 



fli - 2a$i /<t>i + aii c i.i 
fli = alici.j 


^61^63 = 364/4*3 + ^ 65 ^ 2^3 + 361 (364 e l,4 +a 65 e l,5) 
^51^63 = a 5l( a 54 e l,4 + a 55 e l,s) 


for (m = 5,6) 


fm3 — a m4®4,4 ^ 2a m 43 m 564 5 + 3m5®5,5 


^61^22 = a 22/4*2 + a 23/(^-l 4*2) + 361/4*1 + 3*1 ( a 22 c l,2 + a23el3) 

^ 63^22 = 364/4*3 + ^5/(^2 $3) + a 64( a 22 e 2,4 + a 23 e 3,4) + 365( a 22 e 2,5 + a 23 c 3,5) 

^61 ^24 = a 26/4>4 + a 27/(^34*4) + 32g/(X.44>4) + 36l( a 26ei 6 + 82761 7 + a 2 8ei g) 

^61^25 = 329/4*5 + 32io/(^54>5) + a 2ll/(X«4*5) + 36l( a 29 e l > 9 + a 210 e l,10 + a 2lieiii) 


for (m = 6,n = 3), (m = 5,n = 4) 

(ml ^m2 — 3ml (3n2®l,2 4" 3113^1,3) 

fm3fn2 = 3m6(3n2®2,4 4" 311363 4) 

fml fn4 = 3m6(3n4Ci 6 + an7C|,7 + a^e^g) 

fml fn5 = 3ni6(3n9Cl v 9 4" 3,1106110 + Snl 1 1 ) 


for (m = 6,n = 2,3) , (m = 5,n = 4) 

fm3fn4 = 3^,4 (^646 + ^7647 + a^g) + a^a^fi + a n7 C 5> 7 + a^g) 


where the e; j ’s are given by 
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and 


e ij = 


1 

tL + qj 


qi = <h q2 = <t>2 qa = 

q4=<t>3 05 = <t>3^2 q6=<t>4 

q 7 = X -3 <j>4 q$ = X 4 <j >4 q 9 = <t>5 
qi0 = ^5<l>5 qil =^6<t»5 


ir»* -* 
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The through-thickness integration terms are given by: 

g 22 = R 22 [ B 3 t 3 + 3B3B4t 2 + 3B 4 t ] /3. 

g32 = R 33 /I 26 O 5B 3 t 7 + 35B 3 B 4 t 6 + 84B 3 B 6 t 5 + 63B^t 5 

+ ( 105B 3 B 7 + 315B 4 B 6 ) t 4 +420 ( B 4 B 7 + b£ ) t 3 + 1260 ( B 6 B 7 t 2 + B?t ) 

g6i = R<S6/3 B 2 t 3 + 3B!B2t 2 + 3B 2 tJ 

g 22 g 32 = R 23/30 B 3 t 5 +5B 3 B 4 t 4 + (10B 3 B 6 + 5B^)t 3 + 15(B 3 B 7 + B 4 B 6 )t 2 
+ 30B4B7t 

g 22 g 6 i = RW 6 2B]B 3 t 3 + 3 ( BjB 4 + B 2 B 3 ) t 2 + 6 B 2 B 4 tj 

g32g6i =R36/120 4B]B 3 t 5 + ( 15BiB 4 + 5 B 2 B 3 )t 4 

+ ( 40B]B 6 + 20B 2 B 4 ) t 3 + 60 ( B,B 7 + B 2 B 6 ) t 2 + 120B 2 B 7 t ) 

g5ig42 = R 45 /I 2 O ^ 6 B(B 3 t 5 + 15 ( BiB 4 + B 2 B 3 ) t 4 +20 ( BjB 6 + B 3 B 5 + 2B2B 4 ) t 3 
+ 60 ( B 2 B 6 + B 4 B 5 ) t 2 + 120B 5 B 6 tj 

g 4 2 = RW 6 o[ 3B 2 t 5 + 15B 3 B 4 t 4 + 20 ( B 3 B 6 + B 4 ) t 3 + 60B 4 B 6 t 2 + 60BltJ 

g 5 i = R 55 / 6 O 3B 2 t 5 + 15BjB 2 t 4 + 20 ( B, B s + Bl ) t 3 + 60 ( B 2 B 5 t 2 + B?t ) 

g 22 = S 12 /S 11 ( 6B 4 CSIG1 t + 3CSIG2 1 2 + 2CSIG3 t 3 ) /6 
g 61 =S 16 /S„ (6B 2 CSIG1 t + 3CSIG4 1 2 + 2CSIG5 t 3 ) /6 

622824 = ^22 [ B3t5v(l)£i + B 4 ( DIFNU ) ] 
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g24 = R22 4 / 1 ( QUANU ) 

g 32 g 34 = R 33 / 2520 2B 3 t 5 ( DIFNU3 ) + 21B 4 t 4 DIFNU2 + 42B 6 t 3 DIFNU4 
+ 210B 7 t 2 (DIFNU) j 

g34 = R33t 3 /210(QUANU2) 

g6ig63 = R66 / 3 [ Bjt ( DIFMI2 ) + 3B 2 ( DIFMI1 ) ] 

g63 = 4R66 / 3t ( QUAMI1 ) 

g 6 ig 34 = -R 23 / 60 [ B 3 t 3 ( DIFNU4 ) + 5B 4 t 2 ( DIFNU ) ] 

g24g32 = R 23 / 60 B 3 t 3 ( SUMNU3 ) + 5B 4 t 2 ( SUMNU4 ) 

+ 60B 6 t8v(l)e, + 6 OB 7 ( DIFNU ) j 

g24g34=-R 2 3t/15(QUANU3) 
g 22 g 63 = R 26 / 3 [ B 3 t ( DIFMI2 ) + 3B 4 ( DIFMI1 ) ] 
g 22 g 61 = ^26 [ BitSv(l)e, + B 2 ( DIFNU ) ] 
g 24 g 63 - 2R26 / 1 [ 5v(l)Sn(l)e? + 6v(2)5Ti(2)eg ] 

g32g63 = R 36 / 60 B 3 t 3 ( DIFMI4 ) + 5B 4 t 2 ( DIFMI3 ) + 20B 6 t ( DIFMI2 ) + 60B 7 ( DIFMI1 ) 

g 34 g 6 i = -R 36 / 60 [ B,t 3 ( DIFNU4 ) + 5B 2 t 2 ( DIFNU ) ] 
g 34 g 63 = -R 36 t / 30 [ 5tl(l)e, ( DIFNU4 ) - Srj(2)et, ( DIFNU5 ) ] 
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g 42 g 53 = R 45 / 60 [ B 3 t 3 ( SUMMI3 ) + 5B 4 t 2 ( SUMMI2 ) + 20B 6 t ( SUMM1 1 ) ] 

g 44 g 5 i = R 45 / 60 [ B 1 1 3 ( DIFNU4 ) + 5B 2 t 2 ( DIFNU ) ] 

g 44 g 5 3 = R 45 / 30 [ 5tl(l)et ( DIFNU4 ) - 5r l (2)e t> ( DIFNU5 ) ] 

g 42 g 44 = R 44 / 60 [ B 3 t 3 ( DIFNU4 ) + 5B 4 t 2 ( DIFNU ) ] 

g 42 = R44l / 15 ( QUANU3 ) 

g 5 ig 53 = R 55 / 60 [ Bit 3 ( SUMMI3 ) + 5B 2 t 2 ( SUMMI2 ) + 20B 5 t ( SUMMI1 ) ] 
g53 = R 55 t / 15 ( QUAMI2 ) 

g 24 =Si 2 /S„ l8v(De,t(B 9 Sii +B 3 S, 2 +B,S 16 ) + (DIFNU) (CSIG1) ] 

g53 =S, 6 / (3S U ) [ t ( DIFMI2 ) (B 9 Sn+B3S, 2 + B,S, 6 ) + 3(DIFMI1) (CSIG1) ] 

g 22 g 2 5 = S 22 [ — B 3 DIFNUB ] 

g24g25 = S 22 ( - 6 / 1 2 ( DIFNUB SUMNU ) ] 

g25=S 22 [ 12/t 3 (DIFNU) 2 ] 

g32g35 = S 33 h / 840 [ 5t 3 B 3 SUMNU5 + 28B 4 t 2 ( SUMNU6 ) + 42B 6 t ( SUMNU7 ) 

+ 420B 7 ( SUMNU ) J 

g 34 g 35 = - llS 33 t 2 / 2 10 [ (DIFNU) (SUMNU) ] 


g35 = S 33 t / 35 [ QUANU ] 



822835 = S 23 t / 20 [ B 3 t ( SUMNU7 ) + 10B 4 ( SUMNU ) ] 

832825 = — S 23 / 20 ( DIFNU ) ( 3B 3 t 2 + 10B 4 t + 20Be ) 

824835 = 1 1S 23 / 10 ( DIFNU ) ( SUMNU ) 

832825 = S 23 / 10 ( DIFNU ) ( SUMNU ) 

825835 = -6S 23 / ( 5t ) (DIFNU) 2 
861 825 = - S 26 Bi ( DIFNU ) 

861 825 ~ 2S 26 / 1 2 ( DIFNU ) ( SUMMI 1 ) 

861 835 = S 36 t / 20 ( B 1 1 ( SUMNU7 ) +10B 2 ( SUMNU) ] 

863 835 = - s 36 / 10 [ 1*811(2) ( SUMNU 8 ) - etfnO) ( SUMNU7 ) 
851845 = S 45 / 20 [ (DIFNU) ( 3blh2 + 10B 2 t + 20B 5 ) ] 

853845 = 3S 45 / 10 [ (DIFNU) (SUMMU) ] 

842845 = S44 / 20 [ ( DIFNU ) ( 3B 3 t 2 + 10B 4 t + 20B 6 ) ] 

845844 = -S44/ 10 [ (DIFNU) (SUMNU) ] 
g45 = 6S 44/ ( 5t ) [ (DIFNU) 2 ] 

where 


SUMNU = [ 6v(l)e, + 5v(2)£b ] 



SUMNUB = [ 8v(l) n e, + 8v(2) n e b ] 

SUMNU2 = [ 8v(l)e, + 28v(2)e b ] 

SUMNU3 = [ 78v(l)Et + 28v(2)eb ] 

SUMNU4 = [ 58v(l)e, + 8v(2)£b ] 

SUMNU5 = [ 68v(l) n e t + 8v(2) n e b ] 

SUMNU6 = [ 48v(l) n e t + 8v(2) n e b ] 

SUMNU7 = [ 78v(l) n e, + 38v(2) n e b ] 

SUMNU8 = [ 38v(l) n e, +78v(2) n e b ] 

DIFNU = [ 8v(l)e, - 6v(2)e b ] 

DIFNUB = [ 8v(l)„e, - 8v(2) n e b ] 

DIFNU2 = [ 28v(l)e, - 8v(2)e b ] 

DIFNU3 = [ 5Sv(l)et - 28v(2)eh ] 

DIFNU4 = [ 35v(l)e, - 38v(2)eb ] 

DIFNU5 = [ 2Sv(l)e t - 38v(2)e b ] 

QUANU = [ 8v 2 (l)e t 2 + 8v(l)8v(2)e b e t + 8v 2 (2)e£ ] 
QUANU2 = [ 28v 2 (1)e? - 38v(l)8v(2)e b e, + 28v 2 (2)eS ] 



QUANU3 = [ 28v 2 (l)e? - 8v(l)5v(2)e b e t + 28v 2 (2)eg ] 
QUANU4 = [ 138v 2 (l) n e? +98v(l) n 8v(2) n e b e, + 138v 2 (2) n e£ ] 
DIFMI1 = [ 811 ( 1 ) 6 , - Sr|(2)e t) ] 

DIFMI2 = [ 2811(1)6, - 8tl(2)eb ] 

DIFMI3 = [ 3S-n(l)e, - 8ti(2)et, ] 

DIFMI4 = [ 4871(1)6, - 8ri(2)e b ] 

SUMMI1 = [ 8r|(l)e, + 8r|(2)e b ] 

SUMMI2 = ( 38n(l)6, + 8n(2)eb ] 

SUMMI3 = [ 68n(l)e, + 8n(2)6b ] 

QUAMI1 = ( 8ti 2 (l)e? - e b e,8ri(l)8ri(2) + eg8n 2 (2) ] 
QUAMI2= [ 38n 2 (l)e? + 8ti(l)8n(2)e b e, + 38ti 2 (2)eg ] 




Appendix B: Two-term Solution 
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The expressions for the two term solution are presented below. In these expressions G is the 
shear modulus of the resin layer and is taken to be 0.25 msi. 




l-e^ y 


] bP’z + bS 0 


(B.l) 


- — ( Xe _X4y - e~* y ) 

<t> 


2z 


[ Sr| 12 . 1 (k,l>eJ k) + S-q 121 (k.2)eg t) J - -^-8ni 2> i(k,2)eg° 


o® = <t> ' 


U«l 2 

B?>z 2 /2 + B^z+B^j 
G[xe-^ y - e~* y ] [Sn 12 .,(k,l)e«z 2 /(t«) 2 + Sn 12 , 1 (k,2)eg t) (l - z/(t° c >) 2 ) 


(B.2) 



1 - ^ A-^y _ _L P -^y 

1 X-l 6 x e 

- V. J m 

[ B?>z + B^j 

(B.3) 

-G«J> 

X e '^ y " ^ (1 - X)i to + 1)e ^ y ■ ^-♦ y ][5v 12 (k,l)e{ k) [6z/(t ck) ) 2 - 2/t® j 



+ 5v 12 (k,2)eg t) [6z/(t®) 2 - 4/t®]] 


ag = <t>^Y [e - ** - e _x * y j [b?>z 2 /2 + B$°z + Bg°] 

(B.4) 

+ G 

e H,y + |-(1 - X)ye^ y - e - ** 

K 

y [8v 12 (k, 1)6® [3z 2 /(t°°) 2 - 2z/t (k >] 



+ 8v 12 (k,2)eg° ^3z 2 /(t®) 2 - 4z/t°° + lj 


og « <J> 2 rXe- X4,y - e-* y l ["b5‘ ) z 3 /6 + b£>z 2 /2 + Bg°z B^ 

A, 1 . 


-G 


+ |-(1 - X)( e~* y - y<t>el) + h^c~ x * y 


8v 12 (k, l)e® JzVd®) 2 - z 2 /t°°] 


(B.5) 


+ 8v 12 (k,2)eS t) [z 3 /(t®) 2 - 2z 2 /t (k> + z 



Appendix C: Cross-Ply and Angle-Ply Laminates 



Uniform Extension 
e n =0.1% 


[0 4 /90 4 ]j 

[90 4 /0 4 ] s 

Ply 

Oil 

(ksi) 

O 22 
(ksi) 

0\ 2 

(ksi) 

Ply 

Oil 

(ksi) 

0 2 2 

(ksi) 

Ol2 

(ksi) 

0° 

192.8 

3.196 


90° 

192.8 

3.196 

0.000 

90° 

15.54 

-3.196 

0.000 

0° 

15.54 

-3.196 

0.000 


Uniform Bending 

Kji = 0.1 


[0 4 /90 4 ] s 

[90 4 /0 4 ] s 

Ply 

On 

(ksi) 

O 22 

(ksi) 

O 12 

(ksi) 

1 

on 

(ksi) 

O 22 

(ksi) 

O 12 

(ksi) 

0° 

-57.76 

-0.661 

0.000 

90° 

-4.677 

0.146 

0.000 

90° 

-1.530 

1.542 

RKQj 

0° 

-19.28 

-0.3416 

0.000 



















































































90 ° 


<*23 


(a) [0 4 /90 4 ] s Laminate 



<*23 


(b) [90 4 /0 4 ] s Laminate 


Figure C.2. Contributions to a 33 - Extension Load 













X 3 /h 



<*33 


(a) [0 4 /90 4 ] s Laminate 



<*33 


(b) [90 4 /0 4 ] s Laminate 


Figure C.4. Contributions to O33 Stress at X 2 /b = 0.999 for [O4/9O4], and [90 4 /04] s Laminates 





<*23 


(a) [0 4 /90 4 ] s Laminate 



<*23 


(b) [90 4 /0 4 ] s Laminate 


Figure C.5. Contributions to o 23 Stress at X 2 /b = 0.993 for [0 4 /90 4 ] s and [90 4 /04] s Laminates 






X 2 /b 


(a) [0 4 /90 4 ] s Laminate 



X 2 /b 


(b) [90 4 /0 4 ] s Laminate 

Figure C.6. Contributions to o 33 Stress at 0/90 Interface for [0 4 /90 4 ] s and [90 4 /0 4 ], Laminates 
- Extension 















Appendix D: Quasi-Isotropic Laminates 



All quasi-isotropic 
e„ =0.1% 



Ply 

Oil 

(ksi) 

1 

IK3MI 

Ol2 

(ksi) 

0° 

19.18 

-0.096 

0.000 

90° 

1.455 

-5.368 

0.000 

45° 

4.858 

2.730 

-3.111 

-45° 

4.858 

2.730 

-3.111 




















Solution Parameters For Quasi-Isotropic Laminates 
e n =0.1% 


Constant 

[0/ — 45/90/ — 45] s 

[90/45/0/— 45], 

[45/90/0/-45], 

01 (1/in) 

66.93 

50.02 

54.86 

02 (1/in) 

81.68 

61.75 

59.83 

03 (1/in) 

57.73 

99.81 

47.06 

04 (1/in) 

128.4 

56.23 

70.08 

05 (1/in) 

59.99 

174.5 

95.28 

*1 

5.698 

1.386 

6.725 

^2 

4.223 

3.189 

7.801 

^3 

3.394 

3.661 

12.32 

x A 

8.207 

21.11 

11.90 


6.748 

5.047 

1.490 

x* 

18.68 

0.352 

1.572 

A,X1(T 5 

(psi) 

1.242 

3.703 

6.204 

A 2 X1(T 5 

(psi) 

-4.318 

-3.801 

-1.057 

A 3 X10" 5 

(psi) 

-0.009 

-1.002 

-0.386 

h (in) 

.020 

.020 

.020 
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TABLE D 3 . Solution Parameters for Quasi-Isotropic Laminates - Group 2 


Solution Parameters For Quasi-Isotropic Laminates 

e n =0.1% 

Constant 

[45/90/— 45/0], 

[45/0/90/— 45], 

[45/0/— 45/90], 

<t>i d/in) 

66.66 

mm 

52.32 

<t >2 d/in) 

82.13 

61.39 

47.25 

<l>3 (1/in) 

57.31 

41.75 

67.31 

<!>4 (1/in) 

130.4 

77.69 

43.01 

' <i>5 (1/in) 

72.33 

128.1 

41.81 

h 

5.048 

3.609 

3.227 

^2 

4.178 

7.492 

4.040 

^3 

3.320 

mm 

25.14 

x 4 

7.882 

7.074 

6.388 

^5 

7.223 

1.856 

8.910 

h 6 

12.59 


5.618 

AiXlO -5 

(psi) 

1.234 

0.504 

1.564 

A 2 X10 -5 

(psi) 

-4.843 

-1.834 

-1.935 

A 3 X10" 5 

(psi) 

0.178 

0.978 

-0.025 

It (in) 

.020 

.020 

.020 




















































Solution Parameters For Quasi-Isotropic Laminates 
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TABLE D.5. Solution Parameters for Quasi-Isotropic Laminates - Group 4 


Solution Parameters For Quasi-Isotropic Laminates 
£„ = 0 . 1 % 


Constant 

[±45/0/90], 

[±45/90/0], 

[0/±45/90], 

<h ( 1 /in) 


56.94 

61.75 

<t> 2 (1/in) 

_ 

61.75 

58.81 

<t>3 (1/in) 

_ 

99.81 

97.46 

<t>4 (1/in) 

_ 

56.23 


<t >5 ( 1 /in) 

_ 

90.95 

62.39 

X, 



3.530 

x 2 

_ 

4.203 

4.302 


_ 

5.488 

8.821 

X4 

_ 

15.96 

19.86 

X 5 


1.722 

5.189 


_ 


5.082 

a^kt 5 

(psi) 

- 

1.682 

1.565 

A 2 X10T 5 

(psi) 

- 

-3.639 


A 3 X10 “ 5 

(psi) 

- 

-1.352 

-1.068 

h (in) 

.020 

.020 

.020 











































